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METHODS  OF  DETERMINING'  NATURAL  MODES  AND  FREQUENCIES  OF 


TRANSVERSE  VIBRATIONS  OF  BEAMS 


CHAPTER  I 


INTRODUCTION 


I Much  effort  has  been  spent  in  finding  exact  solutions  of  the 

differential  equations  of  vibration  problems.  In  only  the  simplest  of  struc- 
tures, viz.  the  uniform  beam,  is  it  possible  to  determine  rigorous  solutions 
that  include  the  effects  of  rotary  inertia  and  shear  resistance  of  the  beam. 
For  beams  whose  cross  section  varies  in  a simple  manner  along  the  axis  of  the 
beam  it  is  possible,  though  somewhat  difficult  to  obtain  exact  solutions  that 
include  the  effects  of  rotary  inertia.  For  arbitrary  variation  of  the  mass- 
density  distribution  and  arbitrarily  changing  elastic  properties  along  the 
axis,  even  if  they  are  expressible  in  analytical  form,  it  has  not  been  found 
possible  to  obtain  exact  solutions.  Consequently  the  emphasis  of  this  report 
; lies  in  obtaining  approximate  analytical  methods  rather  than  in  working  out 
the  solutions  of  the  differential  equations  of  a class  of  beams  with,  say, 

I cross  sections  that  vary  according  to  simple  analytical  functions  along  the 

I azis  of  the  beam.  The  differential  equation  of  the  problem  will  be  developed 
e 

| and  it  will  be  seen  that  it  does  not  lend  itself  readily  to  the  determination 
| of  approximate  numerical  solutions.  It  is  much  more  advantageous  to  set  up 
I the  problem  in  its  integral  equation  representation,  because  integral  ecrua- 

jj 

| tions  can  easily  be  solved  by  approximate  methods.  The  reason  lias  in  the 

I 

& 

| entirely  different  approach  that  leads  to  the  integral  equations  in  contrast 

| with  the  aifferential  equation  methods.  The  differential  equation  considers 
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the  olastic  structure  in  its  elements;  it  studies  the  motions  of  a partiela 
of  the  aggregate  that  makes  up  the  structure  and  "builds  up  - "by  integration  - 
the  motions  of  the  structure  as  a whole.  These  muet  be  compatible  with  the 
conditions  of  support,  and  it  is  here  that  the  real  difficulty  and  laborious- 
ness of  the  differential  equation  procedure  enters;  the  solution  of  the  bounda- 
ry-value problem,  or  the  evaluation  of  the  constants  of  Integration,  The 
integral  equation  considers  the  elastic  body  as  a whole  and  it  has,  so  to 
speak,  all  the  boundary  conditions  of  the  problem  written  into  it.  It  is  true, 
that  it  builds  up  the  motions  of  the  body  from  "elements",  but  these  are  not 
those  of  isolated  particles  but  simple  motions  of  the  whole  structure.  The 
integral  equations  are  not  the  only  ones  that  consider  the  body  as  a whole, 
for  the  variational  methods,  better  known  to  the  engineer  as  energy  methods, 
also  formulate  the  problem  for  ths  body  as  an  entity,  and  it  is  convenient  to 
say  that  mechanical  problems  of  vibrations  can  be  stated  in  either  differen- 
tial form  or  in  integral  form.  Here,  attention  is  focused  on  the  integral 
f ormulations . 

A few  comments  have  to  be  made  about  the  term  "exact  solution"  that 
has  been  mentioned  above,  i.e.  it  has  to  be  stated  under  what  set  of  condi- 
tions an  analytical  solution  shall  be  called  exact.  The  equations  of  the 
problem  shall  be  linear.  This  limits  the  magnitude  of  vibrations  to  so-called 
small  vibrations  and  it  requires  linear  stress-strain  relationship.  Further- 
more the  simplified  beam  theory  shall  be  valid.  This  means  that  the  Bernoulli- 
Navier  hypothesis  - plane  sections  which  are  norts.1  to  the  undefcrmad  elastic 
axis  remain  plane  and  normal  to  the  deformed  elastic  axis  - is  valid,  ar.d  thus 
any'  exact  solution  is  not  rigorous  in  the  sense  of  the  mathematical  theory  of 
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CHAPTSS  II 


MATHEMATICAL  PREREQUISITE 


+ a 
LrK. 


(b  C * 
1 ' /« 


(1) 


The  ranges  of  the  indices  i,j,  k,  1 shall  "be  n^,  n^,  n~,  n^,  respectively. 

Matrices  are  designated  by  arabic  capitals;  vectors,  that  is, 
matrices  with  either  one  row  or  one  column  only,  by  lower  case  arable  letters. 
Row-vectors  have  superscripts,  column -vectors  have  subscripts.  Bold  face 
letters  (or  underlined  letters)  shall  be  used  where  confusion  of  matrices 
with  their  elements  is  possible.  The  following  are  equivalent  notations  of 


1.  Hotations.  Matrix  Calculus 

Matrix  calculus  will  be  used  extensively  in  the  following  sections. 
Thus  it  is  necessary  to  list  the  notations  that  are  uaed  in  this  paper.  They 
are  in  accord  with  continental  practice,  offer  the  advantage  of  simplicity  and  | 
economy  of  print,  and  are  closely  related  with  indicial  notation. 

In  indicial  notation  - as  it  is  used  hereafter  - the  convention  is 
established  that,  if  an  index  occurs  more  than  once  in  a product,  a summation 
over  the  index  is  to  be  carried  out  over  a number  of  integers  ranging  from  one 
to  n,  unless  otherwise  specified. 

For  example; 

cub,-  - £«.•*.  ■ «ifv  - E *.■  A>  . a.:ib 

fye  C!K  • S 2?  &v  bji  -tK  = E a.ij  ( Z‘  l . 

a“(6"C'«  ‘ ' A *Kr>.C*4*)*  ■ 
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an  m-row  n-column  matrix,  or  an  a l)y  a matrix! 


I &,p 


.s 


S * (* 


*S  L , 


a 


SZ. 


or  alio  : a*.  = { C,n  Qtr>  ■ " , . 


(2) 


Multiplication  of  a matrix  by  a scalar,  - as  multiplication,  of  every  element 
of  the  matrix  by  the  scalar; 

x.  _A  . ( xo,ij ) (3) 

Addition  and  subtraction  of  matrices,  - as  algebraic  addition  of  corresponding 
elements  of  the  matrices.  (Corresponding  means  situated  at  equal  places, 
having  the  saae  address  or  coordinates) 


A"  B +C  - "D  * (rt/) ) Where  cJ  c j - , 


(4) 


Addition  and  subtraction  are  associative  and  distributive  but  are,  of  course, 
restricted  to  aggregates  of  matrices  which  have  all  the  same  number  of  rows 
and  the  same  number  of  columns. 

Multiplication  of  a vector  by  another  vector,  - the  familiar  scalar  product 
of  vector  algebra  is  written  as; 


{Q,  b : •*  hi  + • 4 £\,  ) “ , 

(5) 
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and  with  it,  the  multiplication  of  matrices  by  vectors  and  matrices  are  easily- 
defined.  The  1 indicates  transpositions  of  rows  and  columns.  (The  dyadic 
product  of  vectors  is  not  needed  here.) 

Multinlication  of  a matrix  by  a column-vector,  - is  a new  column  vector  whos3 
i-th  component  is  the  scalar  product  of  the  i-th  row  vector  of  the  matrix  with 
the  column  vector: 


b ■= 


t 

t 


i wifti  <jrb  „ a,,-  hi 


Multiplication  of  two  matrices,  - is  a new  matrix  whose  ij-th  element  (that 
is,  the  i-th  component  of  the  j-tb  column  vector)  is  defined  as  the  scalar 
product  of  the  i-th  row  vector  of  the  first  factor-matrix  into  the  j-th  column 
vector  of  the  second  factor-matrix.  The  two  matrices  must  be  conrpat ible,  i.e. 
the  number  of  the  elements  in  the  rows  of  the  first  matrix  must  he  equal  to 
that  of  the  elements  in  the  columns  of  the  second.  Symbolically  the  product 
is  defined  as; 


[gt  \ 

a® ' 

1 

* 

• 

4 

1 (kW*’  b„)  = 

i 

\ 

a'b, 

2 d a bv 


= ) 
x 3 / 


l-  ■ , s J 


The  matrix  product  is  associative,  distributive  but  not,  in  general,  cosrauta- 
tive.  The  following  relations  are  true; 


% 


C4 
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diagonal  matrices.  The  leading  (or  main)  diagonal  of  a square  matrix  is  the 
line  which  contains  the  elements  with  tvo  equal  subscripts,  (a^) . It  is 
only  defined  for  square  matrices.  A matrix  whose  elements  are  all  zero  is 
called  the  zero-matrix-  If  at  least  one  of  the  elements  of  the  leading  diago- 
nal of  a square  matrix  does  not  vanish,  while  ail  others  are  zero,  then  it  is 
a diagonal-matrix.  A diagonal  matrix  with  elements  all  equal  to  unity  is 
called  unit-matrix-  A diagonal  matrix  with  elements  all  equal  to  a scalar 
is  called  a scalar -matrix. 

, If)  / 01  •'■■«) 

The  unit  matrix,  \ Q ' ' ' j 

hap  the  property,  A£~  G • (9) 

/«U  0 o\ 

The  diagonal  matrix,  = 0 o | 

\ 0 O • • •'olvJ 

can  either  premultiply  or  poatmu.lt iply  a matrix  ^ , yielding  DA  and  AD 
respectively.  Premultiplication  hy  a diagonal  matrix  multiplies  the  rows  of 
the  second  factor  -,  postmult iplication  by  a diagonal  matrix  multiplies  the 
rows  of  the  first  factor  - by  the  corresponding  diagonal  element,  viz.: 

/cU„  «UW  \ U,a„  d^d.L  d>d,3  ... 

Taa  ... 

4,a„  . Aft  . : 


*>ia,  1 


' d,  a., 

L d 3 4 >3 

ft. Li 

Q'll  ft-]. 3 • * ■ 

i 

■ < «, 

• <h  4. 

Repeated  products  of  matrices  with  the  same  factor  are  written  as 
powers  of  that  matrix.  The  zero-th  power  of  a matrix  is  defined  to  be  the 
unit  matrix: 
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A°= 


:-a*  - a?^ 


Operations  on  matrices  that  are  frequently  used  are  defined  next. 


imposition,  - is  an  interchange  of  rows  and  columns  and  designated  by  a 


prime  ' - Tinas 


The  transpose  of  a product  is  equal  to  the  product 


of  the  transposed  factors,  in  reversed  order, 


* A % 


for,  if  C*  AB  - (aiK  ) --(Cij)  ; then 

C -(ajf-c^i)  . Uj<  b,o  --  (boOK,)'  - B’A*  . 


The  inverse  (or  reciprocal)  of  a matrix,  - is  defined  as  that  matrix 


whose  premultlplication,  or  postmultiplication,  with  A gives  the  unit 
matrix  £ , viz.: 


A“-AA"-E  . 


The  theory  of  determinants  shows  that  the  ij-th  element  of  the  inverse  matrix 


is  equal  to 


(Q'Jl) 


C.f. 

! A l 


where  c.f  is  the 


of  tho  element  > and  det  A^llAl} 


determinant  of  the  matrix 


More  details  may  he  found  in  any  toxt-hook  on  linear  algebra*  An  easy  oxposi-  | 


tion  is  V.  Lt  Ferrarfs,  Algebra,  Claredon  Frees  1941*  and  a very  clear  and 


A rr»A  S"M  k C *T*T  ^ & arr\  -vr-i  * t ^r*ir  s-\  r>  r->  >. 
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more  advanced  treat  tee  is  iormuhl^,  Mat  risen  t Springer  1950,  (in  German) 
The  inverse  of  a product  is  equal  to  the  product  of  the  inverse  of 
of  the  factors  in  reversed  order* 


then- only:  B)  ® A 


Frequently  multiple  mtriz  products  will  be  encountered  which 
originate  from  multiple  summat ions , and  the  problem  of  writing  such  summations 
as  matrix  products  arises.  Thus  a few  typical  matrix  products  will  be  written 
in  summation  form.  Only  the  ij-th  element  is  listed.  The  ranges  of  i and  J 


are  one  to  a and  one  to  n,  respectively. 
1.  Example:  the  simple  summation 


Sij  - on*  bjK  * Zr  <*i*  isfiv&n. 


I?  the  order  of  the  indices  in  the  second  factor  were  to  be  reversed  this 
summation  would  coincide  with  the  definition  of  the  matrix  product  AB  . 

As  a simultaneous  reversal  of  the  indices  and  transposition  leave  the  second 
matrix  unchanged,  (bK^);  , the  following  is  true: 


1 5*  (SO  = «*  (a*  (6*/)’)  ■ AB  . 1 

1 (16)  i 

I If  the  Individual  products  is  the  summation  2d  bKj  I 

1 are  aiultiplled,  or  weighted,  by  factors,  d^  , that  is,  if  the  summation  1 

1 S-;  * Z* Q,h  t>„- of  Is  given,  observe  that  the  weighting  factor  (or  numerical  1 

g £37  w ffj  H,  S 

1 integration  factor)  carrieB  the  column  index  k of  the  factor  . This  f 

I 1 

c\  P 

| means  that  every  column  of  the  matrix  A iQ  multiplied  by  a certain  factor  I 
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| . The  second  of  the  equations  (10)  on  page  7 shows  that  such  a matrix  g 

§ A ~ /<•  o'!  I 

| is  the  product  of  a matrix  A , postmult iplied  by  the  diagonal  matrix  D * ( o 3^  )* 5 


Thus  one  obtains: 


" df  ) = f( (li*  (/^  (b/rj.)  ')  = A 


which  is  also  equal  to,  S = A ( B D)*  ; sines  X)  = D 


2.  Example:  given  the  matrix  S whose  ij-th  element  is  defined  by  the  double 
sum:  S.j  s . The  object  is  to  write  3 as  a matrix  product. 

If  one  sets  piV  - a.-k  hhK  , the  problem  is  reduced  to  that  of  the  first  ex- 
ample, for  then,  S:j  = ^ ^ --  piK  . 

This  is  the  Lj-th  element  of  the  product  PC  , where  P=  A 0 , and,  since  the 
matrix  product  is  associative,  there  follows: 


S * '(a»/bJK)  = /4  8C  . 


If  the  two  summations,  over  h and  k,  are  weighted  with  two  sets  of 
factors,  say,  /^/zir>Kx)  and  A7=  , a double  summation  of  the  following  type 

is  obtained,  £ty  - dt  A C r).  Yf)  , and  it  can  (using  the  result  of  the  pre- 

vious  example)  be  written  as  the  ij-th  eluent  of  the  matrix  product  listed 


below: 


A)  • 


0,  * 0 


| If  tne  order  of  the  indices  does  not  agree  with  the  one  listed 

I above,  transpositions  must  be  made.  Not  ail  summations  can,  of  course,  be 

I written  as  tvnical  elements  of  matrix  nroducts:  the  indices  of  the  factors 

P 

J & r*i  a tr«  a a. trv  m n?  tmrr  mvr  nror  & n/^i  x 

? 2 i*w  oir\*jL:uiu,o 
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must  be  such  that  they  can  be  arranged  so  that  the  first  index  of  a factor 
is  the  same  as  the  last  index  of  the  proceeding  factor.  The  number  of  pcssi 
ble  weighting  factors  is  the  same  as  the  number  of  summations.  So,  for  in- 
stance,  the  quintuple  sum:  Su  - Q,  bcdfif  'trtopd  , where 

mk’  nl*  °m’  Pa*  0X0  weighting  factors,  can  be  written  as  the  ij-th  ele- 

ment of  the  matrix  product:  ^ ® ( S;i)  * AMBNCODEaF. 

Summat ions , as  they  have  been  illustrated  in  the  two  examples,  will  be 
encountered  when  certain  integrals  are  evaluated  numerically. 


Another  concept  of  greatest  importance  and  frequent  application  is 
the  problem  of  the  series  expansion  of  arbitrary  functions.  It  will  be  re- 
viewed briefly,  in  somewhat  more  general  form  than  it  is  found  in  the  usual 
texts  of  mathematical  physics,  such  as  H.  & B.S.  Jeffreys,  Mathematical 
Physics,  Cambridge  University  PresB  1950,  or  E.  C our ant  & Hilbert,  Methoden 
der  mathematlschan  Physik.  Xnterscience,  New  York. 

Formulation  of  the  problem: 

Given  a set  of  non-overlapping  sub intervals  of  the  line,  i.e. 
intervals  of  the  form:  x In  this  are  defined  a piecewise  contin- 

uous function  f(x)  and  a set  of  functions,  designated  by  ^ ^C*3  ^ . All  the 
functions  are  subject  to  the  following  restrictions: 

(1)  They  must  b©  integrable  square,  i.e.  their  squares  must  be  int ©gr- 
ab Is  in  the  interval. 

(2)  All  inner  products  of  the  functions  must  exist  in  the  interval. 


Bacall  that  the 


of  two  functions  g(x)  and  h(x)  over 


the  interval  a x 6 b is  defined  as: 


i e.t./'v  in  a cti,‘  iwr\  cTmirniDrc  DrcrADru 

l\^  i 1A  U * AUG  1 U'lVAdJ  * 
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[y*)hz*}4t  3 ($,  h)  = (hj§) 


The  inner  product  of  a function  with  itself  1b  called  th©  gora  of  the  function 
and  written  as: 


= ($,})  = 5<£m  ** 


(18) 


introduce  the  following  abbreviations: 

(?»,£*)•  A -•  A } 

T $w  ) the  /i  -component  of  f in  the  s . 

It  is  desired  to  approximate  f(x)  by  a linear  aggregate  of  the  ^'s, 

f « I £*J  s C*]  • ^9) 

in  such  a way  that  the  approximation  is  best  in  the  least  squared  error  sense, 
i.e.  that  the  integral  over  the  square  of  the  error  is  a minimum,  viz. 


Mi*]  --  “ NTf'f)  ‘(f-f)  f’fj  ” 


(20) 


The  problem  consists  now  in  determining  the  expansion  coefficients 
By  in  such  a way  that  M(x)  takes  on  its  minimum  value.  This  implies  that  all 
the  partial  derivatives  of  M with  respect  to  all  aa  must  be  zero: 


9M  & 0 
/ ^ 


for  all  /h  s 1, 2, 


(20') 


This  set  of  n linear  non-homogeneous  equations  determines  the  • 
Carrying  out  the  partial  differentiation: 

I r 

| and  noticing  that  f'  “ fy*  , there  follows: 


ACDA  C*5  A UXir'  A AJF\  CTDf  ir*rri  inro  nror  i 

nnL/  D ! I\UA  t 'JrVGD  rv  LJLA^1,  -f* 


ranrerat: 


uass  a fi^orBEiantmH 


CONTRACT  NO.  H5  ori  - 07833 

ft e ~2  jtf f*  & --  o . 


Cancelling  the  factor  -2  and  performing 


i the  multiplication  in  the  integrand*  one  obtains: 

ttfr**  - g(k~Q  • (&s  sk  is  a constant  factor) 

| which  can  also  he  written  more  conveniently  as  an  equality  of  inner  products: 


(fif/j)  “ &-'(?»,  fr)  • (20") 

With  the  abbreviations  (18)  this  equation  can  be  written  as: 

€ » ^3  (2U 

where  is  the  matriz  of  the  £>  , and  C and 3 are  the  column  vectors: 

yr* 

c*  {<..a - > a «,j  • 

Equation  (21)  has  the  solution;  ^ * S^*C  , (22) 

Th©  lj«th  element  of  the  inverse  matrix  is  given  by:  (b^)  c f. 

By  virtue  of  the  imposed  restrictions,  c exists  (i.e.  is  finite), 
and  thus  a exists  if, and  only  if,  exists.  This  implies  that  the 

determinant  of  the  matrix  £>  must  not  vanish,  viz. 

11 B j|  ^ o , (23) 

If  one  writes  this  in  expanded  form,  and  introduces  the  definition  of  b^ , 
one  arrives  at  the  fundamental  equation,  the  so-called  Sr&aian, : 


<fn) 


I (ft  <f.) 


ft.  <f< 


If*,  ?4 


^ 0 


which  tne  functions  ^ must  satisfy,  in  order  that  the  expansion  (19)  is 


possible. 


\ f ^ s~*3  a i Mr\  tf'T'i  mrc  nccr  a oru 

i £ V l 1 \J  £SXn*  nLrJljiKUl 


fflvransnnfflcni 
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The  non-vanishing  of  the  Oraaii&n  implies  that  the  functions 
are  linearly  independent,  for,  if  this  were  not  so,  a relation 


d,  (f,  <■  yU  ♦■•••+■  c/n  cfn  = Q (24) 

with  not  all  coefficients  zero  would  he  true.  To  prove  this,  note  that 

from  equation  (24)  any  of  the  f ' s can  be  expressed  as  a linear  combination 
of  all  others.  After  introducing  this  into  the  equation  (23* ) one  obtains 
a determinant  with  a column  that  is  a linear  combination  of  the  other  columns. 
But  guch  a determinant  is  always  zero,  hence  equation  (24)  cannot  he  satisfied. 

Returning  to  the  expansion  problem,  the  approximation  f(x)  (eq.  19) 
becomes  with  eq.  (22) 

f = a>%  * - f'Q'c  . 

Obviously  M(x) , as  a definite  integral  with  non-negative  integrand, 
is  never  less  than  zero.  From  this  remark  a generalized  form  of  the  Bessel 
inequality  and  the  Parseval  equality  can  easily  be  derived. 

M was  given  as: 


M zvcp,)  ^ 0 

Expanding  and  carrying  out  the  integrations,  with  the  abbreviations  (18), 
this  becomes: 

/V7 ' 2^(f,  f.)  *- (%,  = 

I'ff-  -2a,C„  +■  a.u  a.^  t>^  o, 

and  introducing  au  from  (22)  as  du  * 6 , one  obtains: 


V - > + ^ !=•<?. 
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I Sines  b b is  the  product  of  a matrix  and  its  inverse  and  thus  equal 
* y* 

\ to  the  unit  matrix,  one  obtains  finally 

i 

| the  generalized  Bessel  inequality:  cu  ^ 4-  /\T f . (75) 

As  the  right-hand  side  of  (25)  is  independent  of  the  range  of  the  summations 
of  and  /4-  from  one  to  n,  the  Btsnnations  on  the  left-hand  side  can  be  ex- 
tended from  one  to  infinity. 

Row,  there  always  exists  a number  n > 5 such  that  the  squared 
error  K is  less  than  ar  arbitrarily  small  e , i. e . 


M -f[  i-au(p, jVx 


< £ 


if  n > 5. 


In  the  limit  n-*-  -»  M becomes  zero  and  the  Bessel  inequality  reduces  to  an 


equality, 

the  generalized  Parseval  equality: 


A rf  . 


The  equation  (26)  can  also  be  called  a generalized  completeness  relation. 

An  infinite  set  of  lineraly  Independent  functions  is  complete  if  it  permits 
the  approximation  in  the  mean  (least  squared  error  sense)  of  any  piecewise 
continuous  function  to  any  desired  degree  of  accuracy. 

It  should  be  noted  that  the  ecuation: 


Jw.  If  -a* (fb] V.< 


l.i.rn.  ( j - ft)  - 0 


.1.1. m.  means  limit  in  the  mean) 


does  not  imply  the  equality: 


(py 


ngamc; 


i 
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The  latter  is  only  assured,  if  the  sequence  of  the  converges  uniformly  | 

la  the  subset  of  the  line,  because  only  then  is  it  permissible  to  move  the  | 

integral  sign  past  the  summation  sign.  1 

Tory  important  for  later  applications  are  the  statements:  j 

Zi&saa  1:  A piecewise  continuous  function  is  uniquely  determined  by  its  ex- 

pansion coefficients  in  a complete  system  of  functions;  or  also, 

Lemma  2:  Two  piecewise  continuous  functions  are  identical  if  they  possess  the 

same  expansion  coefficients. 

The  proofs  follow  from  the  fact  that  the  difference  of  two  such 
functions  with  the  same  expansion  coefficients  has  an  expansion  with  all 
coefficients  equal  to  zero.  And  thus,  by  virtue  of  the  generalized  Pars aval 
equality,  its  norm  vanishes.  But  this  implies  that  the  difference  itself  is 
zero,  hence  the  two  functions  must  he  identical.  The  expansion  coefficients 
in  a complete  set  describe  a certain  function  uniquely,  even  if  the  series 
expansion  converges  only  in  the  mean  and  not  in  the  usual  sense. 

It  has  been  shown  so  far  that  any  set  of  linearly  independent 
functions  is  a complete  set.  Among  these  the  orthogonal  (and  normalized) 
sets  occupy  a preferred  position,  because  the  determination  of  the  expansion 
coefficients  of  a function  in  an  orthogonal  system  is  greatly  simplified,  as 
will  be  shown  presently. 

Two  functions  are  said  to  be  orthogonal  if  their  inner  product 
| vanishes,  and  they  are  normalized  if  their  norm  equals  unity.  Thus  the  set 
{ 'y/')  is  orthogonal  and  normalized  if  it  obeys  the  two  equations: 

(29) 


orthogonality  relation:  (yC,  Vl)  = f £ /U  " ^ 

) O *'crr 


normalizing  equation: 


NrlyL,  - j 


(3C)  1 


k vnr\  mi  a Awn  cTonmiDcc  occr  a dpu 
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g The  equations  (29)  and  (30)  can  also  be  united  as 

y&)  = s 

in  which  case  the  set  is  said  to  be  ortho-normal. 


With  equation  (31)  the  matrix 


becomes  a unit  matrix  ana  thus 


equation  (’22)  reduces  to: 


a=  C 


which  determines  the  expansion  coefficients.  They  are  the  components  of  f 
in  the  set  of  the  ■njr's. 

The  Gramian  (23')  is  now  equal  to  unity,  and  thus,  orthogonal  and 
normalized  sets  of  functions  are  certainly  linearly  independent.  (This  re- 
mains valid  if  the  orthogonal  set  is  not  normalized,  for  the  Iranian  is  then 
equal  to  the  product  of  the  norms  of  the  members  of  the  set.)  With  3 a 
unit  matrix  B”'  is  also  a unit  matrix,  and  the  double  sums  in  the  generalized 
forms  of  the  Bessel  inequality  and  Pars aval  equality  become  simple  sums  over 


the  c*s  with  equal  subscripts: 

Bessel  Inequality  fci  t-N-f  y 

l 

r 

Parseval  equality  Ct.  = /V* f J 


(33; 

for  orthogonal  functions. 

(34) 


With  equation  (32)  the  approximation  f(x)  is  given  by: 


f = (fw)yc  . u~  ,jZ> 


| which  converges  in  the  mean. 
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raaar  a mjumuamamjoil 


Let  a bean  with  variable  cross  section  be  defined  as  follows i 


x-x  is  a straight  line  and  N is  a plane  normal  to  it.  Q1  is  the  intersection 
of  x--x  with  K . The  boundary  curve  f^(s)  b -s  subject  to  five 

restrictions: 

(1)  0*  is  the  centroid  of  the  area  enclosed  by 

(2)  y^  and  z1  are  the  centroidal  principal  axes  of  the  area- 
moment  of  inertia  of  the  cross  section; 


a non 


VI  AQTTP  AMH 


' a a 


QTD 7 [PTI  IDE'S 

t V A * 


DTOADrW 
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(3)  the  y^  and  Zj  axes  are  always  parallel  to  the  y and  z axes 
respectively; 

(4)  f2(x)  is  a function  of  x only  and  not  of  time,i.e.  the  cross 
section  does  not  change  its  shape  during  vibration,  (if  the 
croea  section  is  hollow,  rigid  stiffeners  must  he  provided  to 
prevent  change  of  shape.) 

The  fifth  restriction  will  he  formulated  presently. 

Suppose  that  some  forces  p(x,t)  are  acting  on  the  beam  in  such  a 
way  that  they  are  in  the  x-y  plane  and  normal  to  the  x-axia.  If  the  beam 
was  Initially  at  rest  it  will  then  vibrate  in  the  x-y  plane  only  when  the 
time  dependent  forceB  p are  acting. 

In  order  to  study  the  motion  of  the  beam,  isolate  a small  slice 
which  is  cut  out  by  two  neighboring  normal  planes  N and  N',  apply  the  inner 
forces  between  it  and  the  remainder  of  the  beam  as  exterior  forces  and  for- 
mulate the  equilibrium  equations.  These  state  that  the  sum  total  of  all 
forces  and  moments,  including  the  Inertia  forces  and  moments,  of  course,  must 
vanish.  Due  to  the  s implifying  assumptions,  only  two  of  the  six  equilibrium 
equations  (in  Cartesian  coordinates!)  are  not  trivial,  namely,  that  the  sum 
of  the  vertical  forces  have  a zero-resultant  and  that  the  sum  of  the  moments 
about  the  z-axis  have  a zero  resultant  moment  vector.  The  following  figure 
shows  the  isolated  slice  of  length  d x.  The  inner  forces  have  been  combined 
into  a shear  force  Q and  a bending  moment  M.  The  indices  n /B  and  fl  r- « 
refer  to  left  and  right.  If  the  beams* s oscillation  is  in  the  x-y  plane  only, 
then  the  resultant  of  th©  shear-stress  intensities  of  the  cross  section  must 
be  parallel  to  the  y-axis.  It  is  known,  however,  that  this  requirement  can 
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only  then  be  fulfilled  if  the  exterior  force  acts  through  the  ahear=c enter 
of  the  cross  section.  As  the  forces  p have  been  prescribed  as  passing  through 
the  centroid  of  the  cross-section,  it  is  thus  necessary  to  specify  that  the 
shear  center  be  in  the  x~y  plane.  It  is  also  Lecown  that  an  axis  of  symmetry 
of  a cross  section  1b  a locus  line  of  the  shear  center.  Hence  the  restriction 
which  must  be  added  is: 

(5)  the  locus  line  of  the  shear  centers  of  all  cross  sections 
nruet  be  in  the  x-y  plane,  i.e.  in  the  plane  of  action  of 
the  exterior  forces  which  is  also  the  plane  of  vibration. 

The  cross  section  must  vary  gradually,  otherwise  the  simple  beam  theory  and 
the  theory  of  the  shear  center  are  no  longer  applicable. 

(Detailed  discussion  of  shear  centers  of  solid  cross  sections  can  be  found 


in  A.  & L.  Foppl, 


vol . 2,  paragraph  78,  Oldenburg  19WU 


Munich.)  Some  other  dynamic  restrictions  will  be  discussed  in  latter  sections. 


J| 


The  figure  shows  the  slice  d x and 


the  resultants  of  the  normal  - and 


shearing  stresses  on  the  left  and 
right  faces.  The  resultants  of  the 


exterior  and  inertia  forces  are 


designated  by  f9  dx.  and  M &3t. 

C is  the  center  of  gravity  of  the 


slice . 


Neglecting  the  dead  weight  whicp  is  an  effect  of  higher  order,  the  sun  of  the 
vertical  forces  must  be  zero; 
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(Xr  - Qt  * p 4 oc  = o 

and  the  sum  of  the  moments  must  he  zero: 

- Mr  + M + ax (Qr  * $e)/2.  * M & 
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0. 


dividing  both  equations  by  4 x and  letting  & x approach,  zero,  one  obtains; 

a'  v 

- A7  -t-  d *■ 

where  the  definitions; 


p 

- 0, 

M 

*0, 

— s M ; ■ - Q,  f Q.r  - Q - $ ) 

ax  0 Ax-+o  0 


have  been  introduced. 

If  only  the  free  vibrations  of  the  bean  are  of  interest,  the  forc- 
ing function  p(x,t)  is  zero,  and  P and  M become  the  d'Alembert  inertia  forces 
acting  per  unit  length  of  bean  axis. 

Due  to  the  postulated  linearity  of  the  problem  it  is  permissible  to 
consider  bending  and  shearing  distortions  separately.  The  total  deflection  of 
the  bean  is  designated  by  y(x,t),  the  deflection  due  to  pure  bending  by  o<(x,t) 
the  deflection  due  to  pure  shear  by  ft  (x,t). 

Thus  7 - ok  -y  (1) 

The  following  two  differential  equations  express  in  terms  of  the  bending 

moment  M and  bending  stiffness  El,  and  ^ in  terms  of  the  shear  force  Q and 
the  shear  stiffness  GAr.  E is  Young's  modulus,  G-  the  shear  modulus, 

G-  s E/  2 (l +'•>),  u s PoiBsons's  ratio  for  lateral  contraction,  I s I_  the 
moment  of  inertia  of  the  cross  section  shout  the  z-,  axis,  Ar  3 A/k  the  reduced 


cross  section. 


The  value  of  the  reduction,  factor  fe  will  be  derived 
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- £1*  = A7  ; 

6Ar/S*  * £ 


(2) 


The  inertia  force  due  to  the  transverse  motion  of  the  team  is  equal 
to  the  negative  product  cf  mass  density  per  unit  length  of  beam  axis,  uCx]  • 
multiplied  by  the  vertical  acceleration  ^:P--  -j* . Before  writing  down, 
the  rotary  acceleration  term,  which  is  equal  to  the  negative  product  of  the 
mass  moment  of  inertia  per  unit  length  of  beam  axis  (about  the  z^-axis)  times 
the  angular  acceleration,  it  should  be  noted  that  the  bending  distortion 
only  rotates  the  beam  slices.  The  shear  deformation  y3  merely  causes  adja- 
cent slices  to  slide  with  respect  to  each  other.  (This  is  the  simplified 
concept  of  the  8simple°  beam  theory!)  For  small  deflections,  the  angle  of 
rotation  is  equal  to  the  angle  of  the  tangent  to  the  beam  axis,  which  is  in 
turn  approximately  equal  to  the  slope  eC.'  . Consequently  the  rotary  inertia 
forces  per  unit  length  of  beam  axis  are  given  by  M -‘fl  . [f  * rt&£sdensih/~] 
With  these,  the  two  equilibrium  equations  of  the  slice  become; 


a-M*  - y I*' 


& 


/n 


0 


^ 0 


(a) 


(fc) 


With 


- j &t\r  t! i mi  note,  e<  from  (JS  j ; 


0.  - K.  4 fXQ[6 Ar 


0 , 


(c) 


Differentiate  (c)  with  respect  to  x; 

& -K!  - ? (x^j*  (f/6j[(x/Arj'a  * (T/Ar)  a’ ) --  o 
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and  substitute  (b)  into  (d) : 


$y+(f/6)[p/4j,&  + CuiMryi]  0 


Express  M by  j as  follows: 


n - - El  ^ - - EX (A£  - fl") 

--  - ex^s  + ex  {a/Ar)V$ 

A7  - - Ely  * l(l*»]I^ (t/Ar)  '&  *■  /ty/Ar]  J 

differentiate  twice  With  respect  to  x ; 

and  introduce  (f)  into  (e) i 


♦ (Ely ' ) t(he)  (^rj/Arf-l(i^)  [(J/Ar)  ’f  G]  - f(l^)  1 F 
+ f [(l/4r)JQ  7^  /Ar]  / 6 » 0 . (g) 

Rearranging  (g) , integrating  (b)  and  introducing  it  into  (g) , the  following 
integro-differential  aquation  is  obtained: 


(EX >{/*i/Ar)fff/6  - 

= 2(w^ 


From  (3)  all  of  the  simplified  equations  listed  in  the  literature  can  easily 


bo  obtained,  for  instance: 
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Bree  transverse  oscillations  of  a beam,  including  rotary  inertia  effects-.  | 

p 

"but  neglecting  shear  distortion.  The  shear-  resistance  is  than  infinite,  i.e.  £ 
Poisson's  ratio  p s -1,  G - go  , and  equation  (3)  becomes:  j 


(tuT  -§ttir 


= 0 


(4) 


This  equation  is  listed  in  E.  Eohenemser  & W.  Prager,  Dynamlk  der  Stahwerke, 
Springer,  Berlin  1933* 

Pree  OBclllation  of  a -prismatic  bar,  including  rotary  inertia  and  shear 
effects.  Setting  El,  yW  , Ar  constant  and  carrying  out  the  differentiations 
in  (3),  yields: 

- (T  4 2.(j+*V  VAr)  1 f f^Ar  ~ Q-  (5) 

This  equation  is  listed  in  St.  Timoshenko,  Vibration  Problems  in  Engineering, 
van  Nostr&nd,  Kew  York,  1948;  where  it  is  integrated  for  some  simple  types  of 
boundary  conditions. 

Assuming  harmonic  action: 

t 

5 (6) 
and  introducing  (6)  into  the  equation  (3)  yields  the  ordinary  differential 
equation  for  the  eigenvalue  problem.  is  the  characteristic  circular 

frequency  of  the  vibration.  The  aquations  (3) , (4),  and  (5)  become  respec- 
tively : 


(Elf  }n  ^ (IV)1  - Vi  ^ * v r$  /*-) * * tf/M  I*  uAr  r 

I Q/Arf  '£(%/&)  (£/K)'&%  dx 

| This  is  the  differential  equation  of  the  natural  modes,  including  rotary 
| inertia  and  shear  terms. 
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(Eitf)"  i-M  ? (r>T  -/« 


O . 


(AO 


This  is  the  differential  equation  of  the  natural  modes,  neglecting  shear 
effects,  but  with  rotary  inertia  effects  included. 

This  Is  the  differential  equation  of  the  natural  modes  for  a prismatic  bar, 
including  shear  and  rotary  inertia  terms. 

Equations  (AO  and  (50  have  been  solved  in  Hoheneaser-Prsger 1 s 
and  1 imoshenko1  s tests  for  the  case  of  uniform  beams  with 'rectangular  cross 
section.  They  show  that  rotary  inertia  effects  become  great  at  higher  fre- 
quencies, because  then  the  nodal  points  are  close  together  and  the  contribu- 
tion of  the  angular  accelerations  is  of  the  same  order  of  magnitude  as  that 
of  the  transverse  accelerations.  It  seems  possible  to  solve  equation  (AO 
also  for  non-uniform  beams,  say,  for  an  exponential  variation  cf  the  height 
of  the  beam.  The  method  of  variation  of  parameters  will  give  an  analytical 
solution  for  the  equation  under  the  assumption  that  X is  Imown.  The  real 
difficulty  arises  when  the  boundary  conditions  have  to  be  satisfied  in  order 
to  obtain  the  characteristic  equation  for  the  eigenvalues.  An  analytical 
solution  of  equation  (3 1 ) » followed  by  the  boundary  value  problem,  appears 
to  be  much  too  difficult  end  tedious  to  be  attempted  here.  Furthermore , 
beams  encountered  in  engineering  practice  have  elastic  properties  and  mass 
distributions  that  can  not  readily  be  expressed  by  exponential  series  and 
rover  series  and  it  is  concluded  that  the  differential  equation  approach  is 
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j not  the  appropriate  one  for  the  problem.  The  integral  equation  methods  and 
variational  methods  will  yield  practical  solutions.  Before  these  are  treated 
the  derivation  of  the  formula  of  the  reduction  factor  MkM  in  A„  s A/k  is 
given. 

The  picture  below  shows  the  side  - and  front  view  of  a slice 
of  the  beam  acted  upon  by  shearing  forces  Q,  only.  Let  it  be  assumed  that  all 
longitudinal  fibres  undergo  the  same  angular  deformation,  i.e.  the  cross  sec- 
tion. does  not  change  its  shape.  Denote  by  the  shear  deflection  of  an 

element  of  length  dx.  One  obtains  an 
equation  for  from  the  fact  that  the 

exterior  work  done  by  C*  during  deformation 
is  equal  to  the  strain-energy  (interior 
work)  of  the  shearing  stresses  produced 
by  Q,.  Q,  undergoes  a displacement  <7/2  and 
increases  slowly  from  zero  to  its  final 
value  Q,.  Thus  the  potential  energy  of  Q, 
is  equal  to; 

4=  (a) 

The  sir  am  energy  ox  the  shearing  forces  is  given  by: 

As  - fYzejfv'i/y  (b) 

v 


The  volume  integral  j"kV/  is  equal  to  dxJ't'c/A  since  d7  is  infinitely  thin 

v 1 

and  X does  not  change  over  the  distance  dx  . The  shearing  stress  inten- 
sity is  given . by: 

Afpruci  AQTin  a wn  ctb!  wTj  noire  orcr  A oru 

4 *4  4*  4 V_/-  M.  A»  4J  « 1 w t a A » 4J.  m 1 \ V_/  w • VW  a.  V ^ « V 4 . 1 J, 
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r * & 5,  //,  A , 


vhar€  S is  the  static  moment  of  the  shaded  area  ahoi^t  the 

2 *L 

Thus,  one  obtains: 

/is  * (ek/ztyftdA  ~(^K/26)f(&St  /ith)  dA  , [dA  - dy-dz] 

"a  * 

and  since  Q,  and  I are  constant  for  a given  x (integration  over  y and  z)\ 

Z 

As  - (Jx/2j)((2/r,ff(sl/ij'w4  . 

Equating  (a)  and  (e):  A » A , there  f 0 Hows : 

Q o 

Qc//}/z  (Mz^&iu)  'S(SWH  VW 

and  solving  for  a : 


One  now  sets: 


djl>  --  (<2dz/ 6 /» /ij'd A . 

A 

As  (A/r,)f(S>/l>)lSA 


The  factor  k depends  on  the  cross  section  only.  With  it  the  differential 
shear  deflection  is  given  as: 

- &d*/ 6^i  • &d*/6Ar . 

which  verifies  the  equation  (2)  on  page  22. 

Values  of  k for  sone  frequent  cross  sections: 

Rectangle:  k s 6 / 5 » solid  circle:  k 2 10/9i  thin—walled  circle:  k s 2.0. 

For  structural  I "beans  one  say  take  with  sufficient  accuracy  the  factor  k a 
1.0  and  the  cross  section  of  the  wet  only:  A/k  a A^^.  ^ ons  takes  the 
total  area,  k varies  from  2.0  to  2.4  for  very  small  and  very  large  sections. 
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2.  Variational  Methods.  Integral  Approach. 
a ) General  Comaent 

In  the  first  section  of  this  chapter  the  eigenvalue  problem  has 
been  formulated,  by  means  of  a differential  equation  and  a set  of  boundary 
conditions.  This  diffe1. entiai  equation  expresses  an  equilibrium  equation 
that  must  hold  for  any  infinitesimal  slice  of  the  beam.  But  this  is  not  the 
only  way  of  deriving  the  differential  equations.  One  could  have  utilised, 
more  elegantly,  Hamilton1 s principle.  For  the  development  of  the  energy 
methods  this  principle  of  least  action  (as  it  is  called)  is  the  natural0 
approach  to  the  eigenvalue  problem.  A very  short  and  neat  development  will  be 
given  which  is  taken  from  G.  Hamel,  Thooretlsche  Mechaaik.  Springer,  Berlin 
194-9,  and  which  does  not  presuppose  a knowledge  of  langrange  * a equations, 

The  derivation  of  the  eigenvalue  problem  follows  an  article  by  E.  J.  M&hly, 

Die  gen&herto  Berechming  von  Eigenwerien  elastischer  Schwingungen  ani sot roper 
Korper,  that  appeared  in  vol.  24  of  Ergebnlsse  der  axakten  Naturvissenschaf ten, 


Let  dm  be  an  inf initesimal  mass  element  of  a system,  w its  vectorial 
acceleration,  &lt@  an  exterior  (impressed)  force  acting  on  it,  and  £ r a 
virtual  displacement  of  the  particle,  r is  the  radius  vector  describing 
the  position  of  the  particle,  v is  its  velocity  which  is  equal  to  the  rate 
of  change  of  r:  v a d?/dt. 


Then, 


;§' s principle  states  that  the  virtual  work  of  the 
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inertia  forces  is  equal  to  the  virtual  work,  * Ae,  of  the  impressed  forces  or 


j £,W>  • 'M?  e ^ d £t  - S /fg. 


(a) 


JcU*,  /v?^r  = j"  or  <J'r  -ST 


0») 


Introducing  the  central  equation,  (b)  into  Lagrange's  principle  (a),  one 

oUalM:  £ &S*  - £T  « SAt 

An  tat egr&t ion  between  t~,  and  yields  Hamilton's  principle  in  its  general 
form: 


J/vw,  vSr  = ] (ST  t S' Ae)  , 


If  the  impressed  forces  possess  a potential,  l.e.  if  £ Ae  is  a 
total  differential  £ ka  m - eT  U,  then  Hamilton'  s principle  appears  in  a 
convenient  form.  Since  the  S-  process  is  independent  of  the  time  variation, 
$ may  he  taken  outside  the  integral  sign:: 


(2) 


ITrom  it  one  obtains  Lagrange's  central  equation  (as  named  by  Heum)  by  writing 

Fi 

dm  wdr  in  a different  way,  i.e.: 

^ £?  s £r  E vSr  -J^£m  V • 

But  if  the  virtual  displacements  are  possible  d ^ r r ^dr,  and 

^£**v*^£  ~ ^ *■  — jCs&w.tf  which  is  equal  to  the  variation  of 

the  kinetic  energy:  0 1,  Thus  there  follows  the  central  equation: 


Cl)  I 


i !“nHr\  r”?  & rTi.^  & airs  pthj  rrv?-^  rvrrr  <v  npi: 
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I L is  the  .tag  rang  ian  function:  L - T - IJ.  The  integral  over  L is  called  the 

r- 

| action,  whence  the  name  of  the  principle.  The  variation  o r can  he  pre- 
1 scribed  to  be  zero  at  and  tg*  Then  one  obtains  the  usual  fora  of  Hamilton1  s 


principle: 


SW  = sf(r-u) 


dt  = 0 


pti 

and  one  says  that  the  integral  ¥ = J L dt  is  stationary  in  the  actual 
motion  of  the  system- 

If  the  system  is  an  elastic  continuum  which  vibrates  freely  (un- 
influenced by  exterior  forces)  then  the  potential  U is  aq.ua!  to  the  strain 
energy  of  the  body.  Thus  one  must  express  the  kinetic  energy  T and  the  strain 
energy  U in  terms  of  the  deformations  (displacements)  of  the  elastic  system, 
introduce  these  into  Hamilton's  principle,  and  thereby  obtain  an  expression 
which  determines  the  displacements. 

b)  Formulation  of  the  Eigenvalue  Problem 

The  free  vibrations  of  a body  are  defined  as  the  simple  periodic 
oscillations  which  are  not  influenced  by  exterior  forces.  Thus,  the  deflec- 
tion of  the  centroidal  axis,  which  is  representative  of  the  deflections  of  all 
the  particles  of  the  bar  (as  has  been  postulated),  can  be  expressed  as: 

y(x,  t ) s y(x)  sin  ( cot  ♦ IF  ) = f ot  (x)  + ^(x)  "l  sin  ( u)fc  ■*-  V')  > (L) 

1 J 

where  y(x)  was  again  divided  into  two  parts,  <x(x)  and  >3(x) , which  repre- 
sent the  bending  deflection  and  shear  deflection  respectively,  viz.: 


at"  = 


(*'  = 
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The  kinetic  energy  cf  the  baas  is  obtained  in  the  following  way.  | 
Consider  a differential  slice  of  the  beam.  Its  length  is  dx.  The  mass-density! 
per  unit  length  of  beam  axis  is  g ^t(x),  the  mass  moment  of  inertia  of  the  | 
beam  about  the  su-axis  (see  figure  on  page  18)  per  unit  length  of  beam  axis  jj 

is  given  by  T- (x) . The  translatory  kinetic  energy  of  the  slice  is  equal 

• ^ *2, 

to  y^y^&xjZ',  and  the  rotary  kinetic  energy  is  equal  to  dx/2.  ( T is 

given  by  It.).  Thus  the  total  kinetic  energy  of  the  bean  is  expressible 


as ; 


r •iJCjif.rp'*)* 


(6) 


The  strain  energy,  in  terms  of  bending  moment  M and  shear  force  Q,, 
is  given  as  U = j[  {M  Zdx/EI  (aidn/cjAr  . With  equation  (5)  this  can 

be  rewritten  in  terms  of  the  deflections  as  follows: 

U * i J(£ IK)2-  + 6Ar((5')Z}d*  (?) 

The  expressions  for  kinetic  - and  strain-energy  can  be  put  into  a form  that 
shows  better  that  they  are  quadratic  forms  in  the  displacements,  viz.: 

r -**[&.&] 

U = £ B [<*.",  o6"]  * |S[(3'>  j3'] 


(8) 

(9) 


where  the  letters  Y and  R are  used  to  designate  the  kinetic  energy  due  to 
vertical  and  rotary  motion,  and  the  letters  3 and  S designate  the  strain- 
energy  due  to  bending  and  shear.  Furthermore  the  notation  is  similar  to  that 
| used  for  the  inner  product  of  two  functions  (see  page  12).  Comparing  (8) 
and  (9)  with  (6)  and  (?)  one  finds  the  definitions; 
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7**  R[(3',  (9'3 1 • cosl(ajt  +w) 

U - ^(b[<  + Sip',  p']j  * 5"?z «• 


(o)t  <•  7|T)  . 


where 


^ = OJ' 


Hamilton’s  principle  states  that  the  variation  of  W must  vanish  for 


any  time  interval  t^  - If  one  chooses  this  interval  to  he  an  integer 

multiple  of  oscillation  periods,  then  the  integrals  over  cos^(cot  + 7/f)  and 
sin^(  o)t  + Y ) are  equal,  and  one  is  left  with  the  variational  equation? 

j[(!K«*)  - SCP',|»'3)  -^(vc-s.y]  * Rl|3'.|»'l)\  « 0 . 0 

This  problem  is  equivalent  to  that  of  finding  the  extrema  of  the 
so-called  Bayleigh  coefficient  K(y)  which  is  defined  as: 


K(vj)  = 


8*S_ 

v+  R ’ 


and  these  are,  by  virtue  of  ( 13 ) » identical  with  the  characteristic  values 

2 . , 
i = COi  a For,  if  one  seta  out  to  find  the  extrema  of  a quotient,  (u/v)  and 

designates  any  such  extremum  by  -t  = (u/v)extr  ’ one  that  the  variation 

of  (u/v)  must  be  zero  at  every  extremum,  i.e. 

S(u/v)  s 0; 

and  carrying  out  the  variation,  one  obtains: 


S(u/i r)  = (vSu  ~ LL&V')//ir: 


S ince  v ^ 0 , 


£(i L/v)  - (^u.  0 ; or  S (il-£v)  = 0’ 

Thus  £ (lL  - X,v)=0  is  equivalent  to  finding  the  extrema  of  X(y) . 
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There  exist  two  methods  by  which  the  quotient  X(y)  can  ba  minimised* 
(i)  One  :a&kes  X depend  on  a series  of  arbitrary  parameters  a, , » 

. . .a^,  and  determines  these  parameters  in  such  a way  that  the 
function  E(y;  a-p.^.a^)  is  a minimum.  This  leads  to  the  re- 
quirement that  n equations: 

3K  . 0 

3a.  ' ° U6) 

be  satisfied.  The  equations  (lo)  are  a necessary  but  not 
sufficient  requirement  for  the  existence  of  a minimum,  and  it 
can  be  demonstrated  that  they  give  an  upper  bound  for  the  high- 
er eigenvalues.  (See  for  instance  L.  Collatz,  Eigenwertnroblem . 
2.  edition,  pg.  27d  etc.)  This  method  is  usually  associated 
with  the  name  of  RIT2. 

(ii)  One  tries  to  make  the  coefficients  of  the  higher  eigen-fre- 
quencies (c^.c-j,..*)  small  with  respect  to  that  of  the  funda- 
mental frequency  c^.  This  leads  to  the  method  of  successive 
approximations  (Iteration).  After  having  found  the  first 
eigenfrequency  one  will  eliminate  its  component  from  the 
Rayleigh  coefficient  and  then  again  try  t-o  make  c^.Ck,  ... 
small  in  comparison  to  Thus  one  obtains  in  succession  the 

first,  recond modes.  The  great  disadvantage  of  this  pro- 

cedure is  the  fact  that  it  is  very  cumbersome  to  ''sweep11  (or 
to  free)  the  Rayleigh  coefficient  of  the  components  of  the 
modes  that  have  already  been  iterated.  As  the  numerical  labor 
incraases  rapidly  for  the  higher  modes  successive  iteration 

methods  will  not  be  treated  here. 
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c)  flie  Ritz  Procedure 

Let  y(x)  s <x(^)f  (3  (x)  be  expressed  is.  the  fora  of  a aeries 

y(x;a^,a2 » • • • >an)  = SI  ai  ^ " aiV'^x^  “ + P‘^3  ^ 

with  yi(x)  = OCtC'x)  + (3t  (x) 

The  functions  yt(x)  hare  to  be  linearly  independent  and  they  must 
satisfy  the  conditions  of  constraint  (boundary  conditions).  Then  B,S,Y,R 
become  positive  definite  quadratic  forms  in  the  coefficients  a^  (see  also 
Collatz,  Elgenwertnroblsme)  and  are  given  by: 


VL},}]  = V[a,a]«  ; "to  Vlj  *!/[>- ^3 

R[^',(3']=  R[a-a]  = rIJa^£1i  J witt>  r‘j=R^*(3^ 

B [<  a"]  = B[d.  = bij  a t aj  ; with  btJ  = B [< . «■  J 1 

S [(3', (3' 3 r S Ca ' a3 s stj a.dj  ; with  sL;  = S[^,(Jj]  . 


From  (1C)  it  is  obvious  what  v.4  sVCjfi.yj]  means;  it;  is  the  Talue 
of  the  following  definite  integral  (to  be  extended  over  the  beam), 

f , „ (18* ) 

v j = !/[*.&]  s ‘ ^ i e‘c’ 

The  matrices  (vy) , (r^j) , (bjj) » (>y)  are  all  symmetric.  Thus  the 
characteristic  values  are  ail  real.  The  Sayleigh  coefficient  K(  , is  an 
extremum  if  the  equations:  dK/**i  - Q are  satisfied,  or  multiplying  both 

sides  with  (V  + R) , one  obtains  the  more  convenient  form: 
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(y  + R)f“  = |^((B^s)-i(V*R)}  = 2i[(ki*^-ifej‘^k  = 0 ■ ( 

For,  if  one  calls  B +•  S - D and.  V + H = E,  one  Las  K - D/H  which  is  the 
function  whose  extrema  have  to  be  determined,  and  it  is  a function  of  the 


The  system  of  linear  homogeneous  equations  (19)  can  have  a non- 
trivial solution  only  if  the  determinant  of  the  coefficient  matrix  is  zero. 
This  gives  the  so-called  secular  equation  (frequency  equation): 


det  u 


<»y] 


where  u^j  = (b  +•  s)jj 


and  » (v  +■  r)  ^ . 


It  should  he  observed  that  the  coefficients  of  BITZ’s  determinant 
(20)  can  be  obtained  directly  from  the  coordinate  functions  y^(x),  as  is 
evident  from  equation  (18),  and  it  is  not  necessary  to  carry  out  the  differen- 
tiations with  respect  to  the  a^'s.  As  the  matrices  (u^)  &Ed  (fc^j)  are  sym- 
metric and  (tjj)  is  positive  definite,  all  roots  JL  ^ are  real  and  they  can 
be  ordered  according  to  size: 

l*  J,  < ...  4 7,  . 

By  a known,  theorem  of  algebra,  the  ratios  between  the  coefficients 
ar^  that  belong  to  a Hitz  value  are  equal  to  the  ratios  of  the  cofactors 
of  any  row  or  column  of  u.  , - t,i  . 
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The  corresponding  solution  functions  Yy(z)  (or  approximations  to 
the  natural  mode  shapes,  as  they  are  usually  called)  are  given  by; 


U (x) 


Ti 


(21) 


The  determination  of  the  is  equivalent  to  the  simultaneous  transformation 

of  the  matrices  (ujj)  and  (t^)  into  diagonal  form. 

The  following  lemma  is  added  without  proof*  (it  can  he  verified  in 
a ’way  similar  to  that  used  hy  H.  J.  Mahiy  in  vol.  24  of  Er&ebnlsse  der  exakrt  en 
Naturwissenschaf ten) * 

The  Hits  values  are  greater  than  the  corresponding  exact  eigen- 
values A 4,  hut  they  are  not  greater  than  the  next  higher  eigenvalue  + ]/ 
If  one  adds  another  function  y .to  the  Hits  expansion,  which  is  linearly 
independent  of  all  yp . • . ,yn,  then  the  Eitz  values  4 of  this  augmented 
system  can  only  be  smaller  than  the  corresponding  ones  of  the  initial  system 
Ji  hut  they  cannot  he  less  than  the  next  lower  value  A of  the  initial 


system. 


-r 


L 4 


i:  t, 


V-  v >... 

It  is  also  possible  to  prove  here  that  the  characteristic  vibrations 
are  no  longer  orthogonal,  i.  e.  that  the  inner  product  of  any  two  different 
ones,  multiplied  by  the  square  root  of  the  naas  density  distribution  function, 
does  not  vanish  any  more.  This  proof  will  he  given  in  the  next  section,  how- 
ever, where  it  will  be  less  abstract.  The  reader  can  find  criteria  for  the 
convergence  :f  the  procedure  in  the  paper  by  H.  J.  Mahiy  which  has  a], ready 
been  quoted. 
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d)  The  Choice  of  Coordinate  Functions  for  the  Eitz  Procedure 

It  has  been  shown  that  Ritzss  procedure  requires  that  a set  of  a 
linearly  independent  functions  y^(x)  be  chosen  which  must  satisfy  the  condi- 
tions of  constraint  of  the  problem.  The  choice  of  these  functions  will  affect 
the  convergence  and  the  exactitude  of  the  procedure  very  much.  If  one  ban, 
h>y  pure  chance,  selected  the  eigenfunctions  as  coordinate  functions,  then  the 
method  would  be  exact.  The  closer  the  agreement  of  the  assumed  functions  y ^ 
is  with  the  eigenfunctions,  the  better  the  convergence  will  be.  It  is  now 
■very  easy  to  find  sets  of  linearly  independent  functions  - that  are  even 
orthogonal  - which  satisfy  all  boundary  conditions  and  which  are  very  close 
to  the  actual  eigenfunctions,  i. a . the  eigenfunctions  of  a uniform  beam, 
neglecting  shear  and  rotary  inertia  effects,  vith  the  same  conditions  of  con- 
straint a®  the  one  to  be  analyzed.  (K.  Eohenemser  & W.  Prager,  gynamik  der 
Stabwerke , Springer  1933 » have  discussed  this  at  great  length  and  have  shown, 
by  many  examples  of  structures  with  very  uneven  distribution  of  mass  density 
and  moment  of  inertia  per  unit  length,  that  the  so-called  uniform  beam  modes 
are  the  best  coordinate  functions  to  use.) 

Tables  for  uniform  beam  modes  can  be  found  in  many  text  books, 
notably  Lord  Rayleigh.  Theory  of  Sound,  and  in  the  text  mentioned  above.  ’Ey- 
tensive  tables  are  added  to  this  report. 

Equation  (4)  shows  that  it  is  also  necessary  to  know  the  bending 
deflection  and  shearing  deflection  separately,  and  thus  the  problem  arises 
as  to  how  to  obtain  Q6(x)  and  (3  (x) , if  their  sum  y(x)  together  with  the 
elastic  properties  of  the  beam  are  given.  Proa  equation  (5)  ona  notices  that 
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"both  M and  Q have  to  he  obtained,  in  order  that  ol  and  |5  can  be  obtained 

by  integrations.  In  terms  of  a loading  intensity  p(x)  per  unit  length  of  beam 

axis,  Q,  and  M are  given  as  first  and  second  integrals  of  p: 

„%  * -x 


Q,  ~ -Jpfifx  ; M = iQdx  = - dr  jf)4x 


Differentiating  (h)  twice  with  respect  to  x»  together  with  (5)  and  (22),  one 
obtains  the  following  general  relation  between  deflection  y(x)  and  load  inten- 
sity p(x) : 


X * ^ % % - % 

J(x)  ~ “ }(d*/6A,)jpdx 


which  reduces,  for  prismatic  bars,  to 

ji  • fib; 


(23') 


(See  for  instance  S.  Timoshenko,  Strength  of  Materials,  vol,  1,  pg.  171) 

The  solution  of  (23)  for  p(x)  with  a given  y(x)  would  be  a very 
cumbersome  undertaking  and  could,  in  general,  not  be  carried  out  analytically. 
But  remembering  that  the  Ritz  procedure  requires  only  functions  that  satisfy 
the  boundary  conditions,  and  that  the  uniform  beam  modes  are  utilized  because 
they  insure  rapid  convergence,  one  can  modify  the  problem  slightly  and  obtain 
OL  (x)  and  (3  (x)  by  the  two  following  procedures; 

(i)  Take  the  mode-loads  of  the  uniform  beam.  These  are  the  static 
loads  that  produce  the  node  deflections.  From  these  one 
obtains  the  bending  deflection  CC  and  shear  deflection  [3  . 
(The  first  eight  mode-loads  for  uniform  beans  with  various 
boundary  conditions  are  tabulated  in  this  report.) 
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(ti)  If  uniform  beam  coda-load. s are  not  available,  estimate  the 

location  of  the  nodal  points  and  apply  alternate  uniform  loads 
between  these  nodes.  Then,  calculate  the  bonding  - and  shear- 
ing deflections  under  these  loads  and  use  them  in  the  Hits 
expansion. 

The  following  picture  illustrates  how  the  third  of  the  linearly 
independent  functions,  y,(x)  s 0C3(%)  -t  jS^x)  can  be  obtained. 


um 

P\ 


0C(x) 


pe*) 


(The  load- intensities  p^.p^.p^  can  be  equal,  or  can  be  determined  such  that 
the  points  Eg  and  do  not  undergo  any  static  deflection.) 

The  underlying  idea  of  this  procedure  is  that  these  deflection  lines 

will  roughly  represent  the  natural  modes  and  thus  one  is  assured  that  the  EITZ 

D -\ 

values  are  close  to  the  actual  and  that  the  procedure  converges 


rapidly.  The  expansion  coefficients  (or  -part leipat ion  factori 

linearly  independent  function  y (x)  in  the  r~th  natural  mode: 

T 


of  the  r-th 


Tr(x)  = ari  yi  (x) 

will  then  he  such  that  the  factor  is  much,  larger  than  the  other  * s 
(where  j ^ r) . This  means  that  the  matrix  (a^.)  is  almost  a diagonal  matrix. 
It  is,  however,  not  necessary  to  choose  the  n linearly  independent  functions 
in  this  way,  and  one  could,  for  instance,  also  seioct  the  deflection  lines 
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due  to  a point  load  that  is  stationed  successively  in  n different  points. 

(It  is  interesting  to  know  that  the  n-th  natural  frequency  of  vibrating  tease 
is  completely  characterised  by  the  location  of  the  nodal  points.  In  beams  with 
both  ends  supported  one  must  not  count  one  of  these.  The  proof  is  given  by 
K.  Hohenemser  & Pragsr  in  SAW',  1931  > vol.  11,  #2) 

The  RXTZ  procedure  offers  the  advantage  that  it  provides  the  analyst 
with  relatively  good,  approximations  to  the  frequencies  and  the  mode  shapes.  It 
has  the  disadvantage  that  the  error  in  the  approximation  cannot  be  estimated, 
and  also  that,  if  one  adds  another  function  y „ to  the  set  of  linearly  in- 
dependent  function  { y } , all  the  labor  of  solving  the  frequency  equation  (20 ) 
has  to  be  repeated.  It  has  been  mentioned  that  this  addition  of  another  coordi- 
nate function  can  only  improve  the  results. 

3.  Integral  Equation  Methods 
a)  Influence  Coefficients 

Consider  a beam  as  defined  in  section  1 of  this  chapter.  If  this 
beam  is  acted  upon  by  a concentrated  load  P which  is  perpendicular  to  the  beam 
axis  sad  located  at  a point  designated  by  then  it  produces  at  a point  "i” 

a deflection  and  an  angular  rotation  0 ^ . Similarly,  a couple  at  8jn 

causes  a deflection  "^ij  ^ andTQdar  rotation  j at  ni”.  The  concen- 
trated load  and  the  couple  shall  both  have  the  magnitude  one. 

Formulas  for  these  influence  coefficients  can  easily  be  derived 

, (and 

moment  and  shear  force  in  the  beam  due  to  a unit  force  (and  a unit  couple)  at 
nkl(.  Then  th.6  above  listed  influence  coefficients  are  given  by  the  following 
integrals  to  be  extended  over  the  beam: 
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« J(MiMj/£j)^x  + j(QtQj/(5Ar)^x 

0ij  * ^JC^iQj/^eJdx 

l9u  = J dx  + |(Qt  Ql/g  Ar)dx 
T)r4j  a J~(R  4Mj/El)dx  + J(Qt  Qj/6Ar)dx 
They  obey  Maxwell's  reciprocal  laws: 
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\ (1) 


J 


Sij  * &JI:  ViJ  s = *^i  ' 

as  is  evident  from  (l).  If  the  point  8i"  has  the  coordinate  x^  - x and  "j" 
the  coordinate  x,  a £ , and  if  both  x and  f are  considered  to  he  inde- 

t) 

pendent  variahles,  then  the  four  influence  coefficients  become  influence 
functions. 

b)  Influence  Functions  (Green's  Functions) 

Let  the  influence  function  corresponding  to  the  coefficient  5.j  he 
designated  by  g.  It  is  clearly  a function  of  x and  § s g a g(x‘,  £ },  and  in 
view  of  Maxwell's  law  it  is  symmetric  in  the  variables,  hence: 

g(x;  % ) = g(  | ; x)  . 

The  determination  of  g(x;  ^ ) is  sufficient,  for,  the  remaining  influence 
functions  can  be  obtained  by  partial  differentiations. 

Lajma  1:  The  influence  function  of  angular  rotation  due  to  a point  load  is 

equal  to  the  partial  derivative  of  g(x:  ^ ) with  respect  to  x. 

AERO -ELASTIC  AND  STRUCTURES  RESEARCH  g 


PAGE  42 

The  following  sketch  illustrac.es  thl 


CONTRACT  MO'.  ^5  ori  - 07833 

In  order  to  obtain  the  deflection  in- 
fluence due  to  a 'unit  moment,  observe 
that  a unit  moment  can  be  considered 
as  the  limiting  case  of  a couple, 

P Ai  = 7fl  » which  is  subjected  to  the 
restriction:  ta.  (P h$)  = <fft  . 

A f-^O 

The  deflection  at  x due  to  a load  ? at 
| * A | is  given  by  g(x;  ^ ) P» 

and  thus  the  deflection  at  x due  to  the 
couple  P Ag  is  equal  to  pjg(x;  £ f d$) 
- g (x;  $ )|  . With  P this 

takes  the  form:  %>{«(*;$+**) -£(*;*)}/£$  . 
In  the  limit  d £-*■  0 the  deflection  at 
x due  to  a moment  % at  £ i 


».  ^ flCx-.g+Af)-  S(*;£)  _ ^ dg  (x ; 0 

3? m 

consequently: 

Lflmntft  Zi  The  deflection  influence  function  due  to  a unit  moment  at  ^ is 
equal  to  the  partial  derivative  with  retpect  to  £ of  the  deflec- 
tion influence  function  for  a unit  force  at  $ . 


Finally: 

Lemm  3:  The  influence  function  for  the  angular  deflection  due  to  a unit 

moment,  is  sque.1  to  the  partial  derivative  with  respect  to  x and 
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£ of  the  deflection  influence  function  for  a unit  force  at  £ « 
The  following  sketch  illustrates  this. 

It  is  pointed  out  that  the  relations 
just  derived  are  not  limited  to  the  G-reerUs 


functions  for  the  deforest  ions.  They  are  squally -H— 


valid  for  those  of  the  shear  force  and  bending 
moment . 


Lamaa  ki  The  influence  functions  for  the  shear  force  and  the  bending  moment 
of  a beam  which  is  under  the  action  of  a unit  couple  located  at 
x = £■  are  the  partial  derivatives  with  respect  to  § of  the 
corresponding  influence  functions  due  to  a unit  force  at  £ . 

The  proof  follows  from  the  remark  that  one  has  only  to  require  that 
the  limiting  value  of  the  couple  be  finite  and  equal  to  TH>  . 

? Hr.  i 

• 1 The  adjacent  figures  show  first  a beam  with 

a unit  force  at  £ and  the  corresponding 
bending  moment  - and  a hear --diagrams , in  which 
a value  at  x is  narked  M(x;  £ ) and  Q(x;  $ ) . 
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Next  is  the  same  bass  with  a unit  couple  at 
? and  with  the  corresponding  diagrams  for 


. “M'  : I 1 1 11  Mil  M I 


,7,5557!  i1 

i j 1 1 j 1 1 1 1 ; 1 1 j !'l! 


j Hots  that  the  shear  fore©  is  always  the 

j H : ' x-derivative  of  the  banding  moment,  it  is 

l.i.  ! ^ ; -|  ! 

i i ■ , ^ 1 ■ 1 i 1 j 

, thus  only  necessary  to  obtain  the  Green's 

1 ''Hj 

: , function  of  the  bending  moment  due  to  a unit 

force.  The  others  are  then  partial  deriva- 

' ! 1 1 1 1 ! ' 'i  " 1 I ! 1 1 1 l!|  1 1 II  lUiTTTIIillllin  »4„  * 4* 

Jin' W-ff.  11708  °f  “■ 

! 1 1 ■ i | ’■  -lllin  1 1 : ■ 11 

!|l|  I1!  j j:  | ! Sometimes  the  evaluation  of  the  Greenes 

I'  j j j function  due  to  a unit  couple  is  easier  than 

the  calculation  of  the  Greats  function  due 
to  a unit  force i followed  by  partial  differ- 
entiations. Then  one  will,  of  course,  de- 
termine it  directly. 

Formulas  for  the  Green1®  functions  of  uniform  beema  are  compiled  in  B.  .Roark, 
Formulas  for  Stress  and  Strain,  McGrav  Hill,  Hew  York  19^3.  The  following 
example  was  taken  from  this  book. 


Illustration: 


,M,Q,y 

1 — £ 

wrj 


uniform  beam,  simply  supported,  acted  upon 
by  a unit  force,  and  also  by  a unit  couple. 
The  Green's  functions,  in  the  interval  from 
zero  to  ^ , are  given  by: 
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M(x;f)  = x(i-f)/£ 


M 3 - %/A 


One  verifies  here  readily  that  the  quant it Iob  indicated  by  a bar  are  the 
| - derivatives  of  those  without  a bar. 


Suppose  that  the  beam  oscillates  freely.  As  shown  in  section  1, 

every  slice  located  at  x « £ is  acted  upon  by  inertia  loadings  which  are 

composed  of  the  inertia  force  -yu yd$  due  to  the  transverse  acceleration 
y and  the  inertia  moment  - "T&'d§  due  to  the  rotary  acceleration  cc'  • 
Let  the  deflection  influence  function  due  to  a unit  load  at  f be  g(x;  $ ) 
as  has  Just  been  established  in  section  b).  Then  the  load  at  | 

produces  at  x a differential  deflection:  — d £ . 

The  deflection  influence  function  g(x;  | ) is  made  up  of  two  parts, 
the  deflection  due  to  bending  distortion  only;  06  (2;  ) and  the  deflection 

due  to  shearing  distortion  only:  ^3  (x;  | ).  Thais, 

3^;f;  = 0C(xjf)  /3(Xj  %)  . 

According  to  Lemma  «f2  of  the  previous  section,  the  deflection  at 

x due  to  a unit  moment  at  § is  equal  to:  , and  consequently, 

the  inertia  moment  -?ai'd§  at  ? produces  a differential  deflection: 

— q-'x!  at  x.  The  sum  of  these  two  differential  effects  is  equal  to: 
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dy  = - | g(x}f)fLyd$  + Toe' dfj 


One  obtains  the  tot&l  deflection  y(x)  integrating  over  all  the  differential 
effects: 


y(- 


■>« = - 1 


et'djJdC  . 


The  natural  modes  are  the  free  harmonic  vibrations  of  the  beam: 
y(x,t)  ~ y(r)sw(cot  + TjfJ  and,  with  cd1  = ^ . one  arrives  at  the  integral  equa- 


tion for  the  natural  modes  of  the  beam: 

y(x)  - + 3|(^—  ret) 


dot(g) 


(2) 


Equation  (2)  is  a homogeneous,  non-linear  Integra-different lal  equation  of  the 
Fredholm  type.  The  function  (X  (x)  depends,  as  has  been  shown,  in  a compli- 
cated way  on  y(x)  • 

The  equation  can  be  simplified  considerably  if  one  sets:  X{%)  ^ y (x / 
i.a.  if  one  neglects  the  shear  effects  in  the  rotary  inertia  terms.  This  is 
justifiable,  since  for  small  vibrations  and  lower  modes  the  functions  y(x) 

I 

and  0C(x)  differ  only  very  little  and  their  respective  first  derivatives  are  j 

I 

practically  alike.  Furthermore  the  approximation  occurs  in  a term  which  is 
only  a correction  term  and  imuch  smaller  than  the  first  expression  in  the  inte- 
grand which  represents  the  deflection  due  to  transverse  motion  of  the  beam 
only.  (This  question  will  be  discussed  in  more  detail  in  chapter  XY.)  With 
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the  simplication  do./dx  - dy/dx  , equation  (2)  appeal's  as 
fb 

y(x) * ’j  + r ( g ) } A g • (3) 

a 

One  can  eliminate  the  derivative  dy  /d%  in  the  integrand  by  partial  inte- 
gration of  the  second  term  on  the  right  and  obtains  then: 


j /(%)=  ^ j <3(x;f)^(C)^U)^f  " \k[«'V 

.a  «i  J L 


or  also: 


. (3') 


Introduce  the  abbreviations: 

G(r,s)  = 3(*i»^(0-||-[t,(l)?2SJ2]  = J( x-,f))Uf)-  (4) 

H (*;?)  = rfpillii?  , 

df 

and  recall  the  notation  of  the  inner  product  of  two  functions,  then  (3*)  can 
be  written  in  the  compact  fora: 


(6,y)  + H yj| 

L ’ J 


Note  well  that  & is  a function  of  % and  § , a Green's  function  which  is 
?|nt  symmetric,  and  y is  a function  of  I also,  if  it  is  in  a round  bracket. 
This  bracket  is  a short  notation  for  the  definitf  integral  of  the  product  of 
the  factors  inside,  and,  since  G-  is  also  a function  of  x,  which  is  not  affect- 
ed by  the  integration  process,  (G,y)  must  be  a function  of  x.  Similarly  H 
is  a function  of  x and  f , and  y is  again  a function  of  f , whenever  it 


is  multiplied  with  another  function.  The  vertical  bar  1 indicates  that  the 
difference  of  the  values  of  the  function  Ey  for  $ ^ b , and  § = a is 
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to  be  taken.  Sine©  s is  not  involved  in  this  process  of  subtract- Lng  the  lover 
limit  from  the  upper  one,  £7  i is  indeed  a function  of  x,  Consequently, 
both  sides  of  equation  (5)  are  functions  of  x alone.  If  the  simplification 
oC  ^ y is  not  made , equation  (5)  reads: 

1 ' C ^ C*  r -if ; | 

y(*)  = ^]\  3(x,?)/*(f)yW  d*  - \ + > . (5,) 

*“4  a ^ ^ J 

From  equation  (5)  one  proves  now  without  difficulty  that  the  charac- 
teristic functions  are  not  orthogonal,  for,  let  y^  and  y j be  two  eigenfunctions 
and  and  the  corresponding  eigenvalues,  then  the  two  following  equa- 

tions must  be  satisfied; 

yi/hi  = (6>^)  + Hytj  , 

]:/*>  = + H>s|  ‘ 

Take  now  the  inner  product  of  the  first  aquation  with  respect  to  y , also 

t) 

the  inner  product  of  the  second  one  with  respect  to  y4 , subtract  the  second 
one  from  the  first  one,  and  note  that  (yi,y^)  = then  one  obtains: 

(*rt)  (*.>*)  3 ((6.y;)>?j)  (Hyi|;>j)  - ((e,^),^)  - (Hy^y,)  . 

If  one  writes  this  out  in  detail,  it  becomes: 

b k o fc  *> 

3 y,(x)dfdx  - j j G(z  O%(Vy,&)  df  dx 

rt  <?■*>  rb 

+ 1 y/A)d-x  --  H(x)f)yj,f)  yt{x)cU  , 

* i U 

Inverting  the  dummy  variables  in  the  second  term,  and  combining  then  these 
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two  integrals  by  factoring  out  y^(  f )yj(x),  one  gets 


h d>  " 


b b 

J S jci* ; f ) - 6 (f ; *)J  yfi)  'Jj  W 

a 

Vr  (•  r 

( H(t;5)kW^W  ~ 


(M 


S*b 


dr 


J y 

(for  Lfj) 

If  the  'beam  were  supported  at  both  ends,  then  y would  be  zero  at  the  ends  and 
consequently  the  second  integral  of  the  right  hand  side  of  (6)  would  vanish. 
Since  G-(x;  § ) ^ Q-(  f ;x) , one  concludes  that  for  th*  s case  the  orthogonality 
cannot  exist.  If  the  beam  is  clamped  at  one  end,  then  the  second  integral  of 
the  right  hand  side  of  (6)  will  not  vanish,  and,  in  general,  orthogonality 
cannot  exist.  In  a similar  way  one  verifies  from  (5’)  that  the  eigenfunctions 
are  not  orthogonal. 

Only  if  rotary  inertia  effects  are  neglected,  can  one  obtain  ortho- 
gonality with  respect  to  a weighting  function  , for  in  this  case, 

equation  (2)  reduces  to: 


e 


(?) 


With  the  new  function  y =?  jjl  y this  becomes: 

cb 

a 

and  if  , yj  and  A;  , A,  are  two  different  eigenfunctions  and  their 

corresponding  eigenvalues  of  the  problem,  there  follows: 
f 


i (tXV) 
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'bluing  it  with  the  first  one,  there  follows: 

,b  b j> 

,x  =f 

' * ' 

In  section  b)  it  has  "been  shown  that  g(x;  § ) = g(  § ;%);  cona8q.uen.tly  the 


/i  i \ f ^ - *b  b 

IT"  v)L2<<x  = i Uj(*ii)- ijMdfd* . 

^ ^ J J l . J 


(8) 


right-hand  side  of  (8)  is  zero,  and  one  finds,  that  for  1 ^ $ 

d x = 0 ; 


I 


b 


with  and  Jj  this  equation  yields  the  generalized  orthogo- 

nality relation  (orthogonality  with  respect  to  the  weighting  function  y.  ): 

f dx  = 0 (9) 

or  one  may  also  state  that  the  functions  '\fJL  and  ^ are  orthogonal 

in  the  conventional  sense.  With  the  terminology  of  chapter  II,  section  2,  it 
is  concluded  that  the  eigenfunctions  of  the  vibrating  beam,  neglecting  rotary 
inertia  effects,  multiplied  by  the  square  root  of  the  mass  density  distribu- 
tion function  form  an  orthogonal  set,  which  can  readily  be  normalized.  For, 

Lf:  ft)®  then: 

iS  * »€ *) = e;>  • (ea  ■ 1 ’ eij*° 1 tj) 


PAGE  50  CONTRACT  N0.H5  ori  - 07833  | 

Multiplying  the  first  equation  by  y^(x),  the  second  by  y^x),  integrating 

over  x between  the  limits  a and  b.  and  subtracting  now  the  second  fr  m the 

first  one  of  the  altered  equations,  yields: 
b b b 

(t.'iJjv.X^r  = j \ ■ 

Interchanging  the  dummy  variables  z and  § in  the  second  product  and  com- 


(10) 
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which  is  the  statement  that  th©  set 


I is  ortho-normal. 


The  integral  equation  (3).  or  (5')>  can  in  general  not  be  solved 
exactly,  and  the  problem  arises  as  to  how  to  find  approximate  solutions  of 
it*  Assume  that  the  exact  solution  can  be  approximated  by  a linear  aggregate 
of  linearly  independent  solutions  y^(x)  “ or  °^(x)  ~ which  satisfy  indivi- 
dually the  boundary  conditions  of  the  problem. 

This  is  again  a RlfZ  procedure.  So  let: 

y i'-x)  ss  y (%)  s Qj  y}  (x)  = <*j [<*,(*)  + /% M J (n) 

and  introduce  it  into  equation  (5),  which  was: 

y s ^ [ (G»y)  + HyJ  ■ . (5) 

This  yields:  ^ « A £(6,  a3  yj)  + H a^J  j = A a,  [(<3,  + H^jj  <12> 

and  one  requires  now  that  the  left  hand  side  coincide  with  the  right  hand 
side  at  as  many  points  Xj  as  there  are  unknown  coefficients  a^,  i.e.  at  n 
points*  (See  section  f)  for  motivation.)  This  gives  a set  of  n linear  homo- 
geneous equations  for  the  coefficients  a.j 

with  the  abbreviations: 


f ' (j-j)  1 UM 


* [a,  a,,  ■ ,a„ | 
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= (G‘j)=  ( lstx‘'5)  ^[s)  il)  > 

= (H‘^=  , 

= 6 + H 

(13) 

one  obtains  the  follow  lag  generalized  matrix  eigenvalue  equation  for  the 
vector  3 : 

ya  = > La  o*> 

Shis  equation  can  he  reduced  to  an  ordinary  eigenvalue  equation  by  premulti- 
plying both  sides  with  the  inverse  of  the  matrix  L—  : L™  > which  exists, 
since  |S  L||  is  not  singular;  thus  with  1 L ■=  £ , E a = a 

l"  ya  = a a 

and  with  t‘y  = B , where  B Is  the  solutioa  of  the  equation! 

Lx  = y B - X , 

one  obtains  the  ordinary  matrix  eigenvalue  equation: 

B 3 - ^ *3 

This  homogeneous  equation  has  a solution  only  If  the  determinant  of  the 

coefficient  matrix  vanishes: 
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(16) 
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det  _y-  ^ LJ  - j|  y - A Lii  = ° - (17) 

| Sq.UB.tioa  (1?)  is  the  generalized  characteristic  equation,  or  the  lambda  equa- 

!tion,  and,  written  out  in  detail,  reads: 

■■■  , ^-xL 


~ ^ I'M  f y^i  ^ y ‘ , ^/nfl  ~ ^ 

It  is  an  algebraic  equation  in  A whose  coefficient  of  A is  equal  to  II  <L 


and  whose  absolute  term  is  jjy  j| . Thus  it  is  an  equation  of  the  n-th  degree 
in  A , if  the  matrix  La  is  not  singular.  With  this  condition  equation  (16) 
exists  and  one  can  obtain  the  character 1st  is  frequencies  from  the  secular 
equation: 


aE  «o, 


b„-> 

b,„ 

bl2~  ?l 

&zrt 

;=0 

’ 08) 

bn, 

brrt  *** 

« 

i 

c 

X? 

that 

(16)  possesses 

non-trivial 

solutions. 

The  matrix  & Is  best  obtained  by  solving  the  set  of  linear  equations  ( 15) • 
It  is  now  necessary  to  examine  hov  the  matrix  L « aq. (13) » is  con- 
structed. Proa  eq.  (4)  one  finds: 


G(^ii) » $(*;$)/■(%)  - — 


T(t)  tlpl> 


where  g(x;  f ) is  the  deflection  influence  function  due  to  a unit  load  at  ? . 

Squat  ion  (3 ‘ ) states  what  integration  has  to  be  carried  out.  The  ij-th  ale- 
sent  of  @ , (J.  , is  the  value  of  the  following  integral  at  a point  with  the 

1 ,3 

coordinate  x s and  with  the  J-th  coordinate  function  y^(x): 
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The  points  are  taken  at  regularly  spaced,  intervals  over  the  "beam  axis.  Let 
the  point  s a a correspond  to  i s 1,  and  x - b to  i ® n;  i.e.,  let  the  length 
(b  --  a)  be  divides  into  8n-!R  equal  intervals.  Prom  (19)  it  is  apparent  that 
every  Gjj  is  the  difference  of  tvo  values,  namely  the  value  of  the  integral  at 
the  upper  limit  minus  the  value  of  the  integral  at  the  lower  limit.  But 
is  a similar  difference,  vizc 


and  consequently  (L 4 and  Ej  < can  be  combined  to  yield  A 4 which  appears, 


written  out  in  full*  as : 

l.  = Aj) ; Aj  = r *T®^r-V)| . (zd 

L . * Jf*a 

The  procedure,  outlined  so  far,  ie  known  as  collocation  using  assumed 

functions.  It  is,  of  course,  not  the  only  possible  method  for  determining  the 

coefficients  Bi  as  well  as  possible. 

J i 

I 

Other  Methods  for  Determining  the  Coefficients  a. : 

3 i 

A more  accurate  result  is  obtained  - but  with  more  work  - if  one 
requires  that  not  the  ordinates  of  the  approximating  function  be  equal  to  the 
values  calculated  from  equation  (12)  at  n discrete  points  (usually  equally 
spaced),  but  that  the  integrals  over  the  n subintervals  (usually  of  equal 
length)  be  equal,  i.e.  the  integral  over  the  k=>th  sub  interval  of  the  left-hand 
side  of  equation  (12)  shall  be  equal  to  the  corresponding  one  of  the  right-hand 
side  of  this  equation.  The  integration  process  has  the  effect  of  making  the 
mean  values  of  the  two  functions  of  the  left-and  right-hand  side  equal  in  n 
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S sub intervale.  Formally  this  yields  the  following  set  of  n linear  homogeneous 

1 

i 

I equations: 


aj j ^ = Aaj]  [(Q*yi)  + H5i| 


dx 


(22) 


Xi- 


l-l 


where:  Z,  “ Xi-i  = 8 -a  , and  usually  - (b  - a)/n  9 

One  can,  of  course,  write  the  values  of  the  integrals  over  the  n 
sub  intervals  as  elements  of  matrices, 


C<j  s 


■r« 


; ^PC 


dtj  m |[C<3^j)  H Vj|]dx  , 


(23) 


*1. 


fc'l 


and  put  (22)  into  the  form:  C a = > Da 

or,  Fa=^a  witt  F = DC  . 


(24) 


The  frequency  equation  associated  with  (22)  or  (24)  is  given  "by: 


IC-^D 


= o 


0? 


F-  a E 


= o 


(25) 


Still  another  method  for  determining  the  factors  is  one  based 
upon  lemma  2 of  chapter  II.  It  was  shown  there  that  a function  is  uniquely 
determined  by  its  expansion  coefficients,  particularly  that  the  difference 
of  two  equal  functions  (called  the  zero-function)  must  have  vanishing  expan- 
sion coefficients.  From  equation  ill-22)  it  follows  that  these  are  zero,  if 
the  components  of  the  zero-function  with  respect  to  all  the  members  of  a com- 
plete set  of  functions  vanish. 

i 


(For, 


‘lC  =0,  or  a„  = h'J’cy.  rf  c^s( f,%)aG,  since  H’ A Q.) 
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In  the  present,  case,  it  is  desired  to  maxe  equal  to  + HyJ  j 

or,  which  is  the  same,  their  difference  equal  to  zero: 

ajYj  ~ ^{(G.y,)  - H^lj  = 0. 

Thus,  if  ( $/*>}  is  a set  of  linearly  independent  functions, 

C -*^{(6^)+  Hyjj  f ^ J 33  0 , i r L2,  1 . . ,Q,  ,1  • ( 2a j 

must  be  satisfied  for  every  i,  and  the  have  to  be  determined  accordingly. 
This  equation  can  also  be  written  as 

<%(», <fi)  = * <*) - Hj5l  , fi)  <«> 

where  i and  j must  have  the  same  range  in  order  that  the  a^’a  can  be  deter- 
mined. If  both  [ y,  ] and  ( ft]  were  complete  seta  (with  infinitely  many  mem- 
bers), this  procedure  would  give  the  exact  solution,  provided  the  series  a^  y^ 
were  absolutely  convergent.  If  only  the  firBt  n terms  of  the  sets  are  taken, 
one  obtains  an  approximation  to  the  exact  value. 


Setting; 


hM  *(»>?;)  > 


= (Ce.yj)  + Hyj  , f>.) 


equation  (2?)  can  be  written  ae: 

= . (29) 

The  particular’  choice  of  the  complete  set  has  a decisive  influence 
upon  the  quality  of  the  approximation.  The  process  Just  described  is  usually 
called  the  method  of  weighting  functions,  and  it  is  said  that  in  eq.  (27)  the 
inner  products  of  both  sides  with  the  weighting  functions  pi  must  be  equal, 
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An  easy  set  to  work  with  are  the  functions, 

x°  , x1,  s2,  . . . . , Xs-1;  (30) 

which  requires  that  the  areas  of  the  curves  represented  by  a.y , and 
}\  ' (G,^)  + H^-jr  with  the  x-axls,  and  the  first  (n-1)  moments  of  these 

areas  shout  the  y-axis  are  equal. 

Another  convenient  set  is  that  of  the  coordinate  functions  used  in 
the  HITS  expansion, 


71’  72‘  y3’  ’ * • * • Tn  J I31J 

and,  if  they  are  orthogonal,  H of  equation  (29)  reduces  to  a diagonal  matrix 


whose  elements  are  the  norms  of  the  y^’s 

/ d,  o o\ 

Dft  - ^ K3  with  D = ( : : **.  , ) > ^Vi»yj)s^S  • (32) 

\0  0 •••«»/ 


Premult iplying  both  sides  hy  the  inverse  of  the  diagonal  matrix;  O 

I 

and  setting  A - ^ , there  follows; 


DKa  = ? a , ^ rvia  = \ a (33) 

where  M - D K . The  inverse  of  a diagonal  matrix  is  simply  the  matrix 
obtained  by  replacing  the  diagonal  terns  by  their  reciprocals,  and  is  there- 
fore again  a diagonal  matrix;  . ( ^d,  " 0 \ 

»'-(i  ■ 

Recalling  the  multiplication  rule  of  diagonal  matriees  with  matrices  (eq.10, 
page  6) , it  is  se9n  that  the  matrix  is  obtained  from  the  matrix  K by 


dividing  every  row  vector  of 


by  the  corresponding  element  of  the  matrix 
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Written  out  la  detail: 


j K/a- 


i K,/d, 


fc'.z/d  i 

&iz/dz 


'■^nydn  &nzj>dn  4 * * ^nn^/otn 

- ij 


(y.  >%)  * -*j 


The  order  of  the  characteristic  values  ( ^ - /’X  ) i-B  aov  reversed, 
One  could,  of  course,  also  have  used  an  ortho-normal  set, 

e< 

la  the  RXT’Z  expansion  and  would  then  have  arrived  directly  at  eq.  33  • 

These  last  methods  are  sometimes  called  methods  of  weighting 

functions* 

The  methods  presented  so  far  are  based  upon  the  assumption  that  the 
integral  equation  (2),  page  46,  (which  is  to  be  solved)  can  be  put  into  the 
form  (5)  in  which  the  first  derivative  of  the  unknown  function  does  not  appear 
(see  eq.  3s) » ana  where  the  linear  i cat  ion  QC  has  been  performed.  The 

elimination  of  V'  was  achieved  by  a partial  integration,  and  in  carrying  this 
out  it  was  necessary  to  calculate  the  partial  derivative  with  respect  to  6 
of  the  Green's  function  H(x;  £ ) , as  defined  in  eq.  4-,  pg.  47=  In  applications 
it  is  often  difficult  to  perform  this  differentiation,  and  in  cases  where  it  is 
easier  to  obtain,  the  influence  functions  numerically  in  the  form  of  a matrix 
(which  is  merely  a function-table  listing  the  deflections  at  a certain  point 
due  to  a unit  load  at  various  other  points  or  the  same  point)  one  would  have 
to  perform  the  differentiations  numerically  also.  These  are  relatively  cum- 
bersome and,  as  with  all  differentiations,  represent  a loss  of  accuracy.  Thus 
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it  is  desirable  to  treat  equation  (2)  directly  without  linearising  it  and 
without  transforming  it  Into  the  form  (5 1 ) < page  4? • This  requires  a modi- 
fication of  the  uethods  that  have  been  outlined,  as  is  shown  in  the  section 
below. 


solutions 


i@  aon= 


The  equation  (2)  of  section  c)  (pg.  46) ; 
where  g(x;  £ )>  Sg(x;  f )/3£  , yi(0  . T(|)  are  known  functions,  has  to 


be  solved. 


Let  the  unknown  function  y(x)  * 04  (x)  + j 3(x)  be  approximated  by 


a EITZ  expansion: 


y (x)  ~ y = *j[otj(x)  ▼ (*)J  > 


y»  = + ft  ? J 

where  the  coordinate  functions  y^Cx)  individually  satisfy  the  boundary  condi- 
tions, are  linearly  independent,  and  are  preferably  rough  approximations  to 
the  natural  modes.  Take  them,  for  instance,  as  equal  to  the  deflections  of 
the  beam  due  to  the  mode-loads  of  a uniform  beam  with  the  .same.,  boundary  condi- 


tiona . 


The  static  mode-loads  of  non-uniform  beams,  neglecting  the  influences 


of  shear  and  rotary  inertia,  satisfy  the  equation  (see  eq.  4*,  chpt.  Ill): 


(Elyff  = ulfty 


K > 


& « 1,2, .... ,n, .. . 
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The  static  deflection  of  the  earns  beam  due  to  a load  p(x)  is  given  by: 


(SIya ) 1 


Comparing  this  equation  with  the  preceding  on©,  it  is  apparent  that  the 
natural  modes  are  identical  with  th©  static  deflections  due  to  the  loads 


Pk  - 

for,  the  equality  of  the  right  hand  sides  of  both  equations  necessitates  the 
equality  of  their  left-hand  sides.  Th®  loads  p^.  are  the  so-called  s^aL^S. 


It  has  been  shown  that  the  function  y%\  are  orthogonal. 

From  the  last  relation  it  follows  therefore  that  the  function  | py j are 
orthogonal.  They  can  readily  be  normalized. 

For  a uniform  beam  with  both  ends  supported  the  modes  are  trigono- 
metric functions  (sine  waves);  other  support  conditions  yield  products  of 
trigonometric  and  hyperbolic  functions.  If  the  uniform  beam  is  gradually 
transformed  into  a non-unifora  beam,  the  modes  transform  (as  do  the  natural 
frequencies)  continuously  into  those  of  the  non-uniforss  bean.  Their  shape 
still  resembles  the  sine  curves  (or  modified  sine  curves)  of  the  uniform 
beam.  The  mode  loads  of  the  non-uni form  beam,  however,  being  products  of 


the  normal  modes  (or  nqa 


as  they  are  also  called)  and  the  mass 


density  distribution  function  fJ-  , will  lock  very  such  different  from  the 

smooth  mode-loads  of  the  uniform  beam.  The  following  picture  shows  on  the 

left  the  (second)  mode  and  the  corresponding  mode-lo&d  of  a uniform  beam. 

On  the  right  is  the  corresponding  mode  and  mode-load  of  a non-uniform  beam 

with  the  same  conditions  of  support. 
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Thns,  it  is  indeed  obvious  that  the  product  of  y,&  times  jj.  (variable)  j 

resembles  yU  yK  very  mu oh,  and  it  is  concluded  that  the  coordinate  func-  | 

tions  in  the  Eitz  expansion  ( 11)  ought  to  be  obtained  as  the  static  deflections 
due  to  the  loads  which  are  the  product  of  the  uniform  beam  modes  times  the 
mass  density  distribution  function  yW  of  the  non-uniform  beam. 

The  evaluation  of  yj  and  °6j  does  not  necessiat©  additional  cal- 
culations, sines  yj  has  to  be  obtained  as  the  sum  of  06j  and 

After  this  excursion  on  a rational  choice  of  the  coordinate  func- 
tions of  the  EITZ  expansion,  there  remains  to  show  how  equation  (2)  is  to  be 
solved. 


Two  basically  different  methods  will  be  discussed  as  in  section  d) , 
vie.  i .collocation  - and  weighting  ^potion  methods.  Since  the  underlying 
principles  have  already  been  established  the  argument  need  not  be  long. 

The  integral  equation  (2),  which  must  be  solved,  is  of  the  form: 


v 

i 


+ (mX)}  , 


(2 * ) 


L 


AERO- ELASTIC  AND  STRUCTURES  RESEARCH 


■*A 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

Department  of  Aeronautical  Engineering 


PAGE  62 


CONTRACT  NO  S5  ori  - 07&33 


therein  the  abbreviations 

vb 


K , y)  ~ f (?)  y (V  4%  , 


w 


I 


have  been  introduced,  and  where  the  definite  integrals  are  written  as  inner 
products.  Introducing  the  two  RITZ  expansions 


" e-'ji 


y * y ~j 

into  (2*)  one  obtains: 


o if  « od  - «,  ot< 


aj7j  ~ ^ { (K>y,)  + (M,  <*j)j 


(11) 


(35) 


The  coefficients  a^  remain  to  be  determined  in  such  a way  as  to  minimize  the 
discrepancy  between  the  left  - and  right  band  side  of  this  approximate  equa- 
lity. Since  a coefficients  &j  are  “free0,  one  must  prescribe  n conditions 
which  have  to  be  satisfied: 

collocation  at  n points  is  one  possible  set  of  conditions. 
orthogonality  to  the  first  n members  of  a complete  set  of 
functions  is  another  one. 

Collocation  at  n distinct  points  transforms  the  approximate  equality  (35) 

Into  the  following  system  of  homogeneous  linear  equations  for  the  a^!a: 


«,•*•(*>  = M ((«[**], y/tf)  1 (m£*:h, «;«];)  ■ 


Al 


/ - s s 1 

\j)Oj  - 
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With  the  matrices  and  vector  defined  by: 


YGyij);  K = K) 


PAGE  63  | 


(%)  ; i = (^}  ; 


(3?) 


'jij  - ; Kcj  ~ kuf3  3 yjCc])  ; mLJ  , \ 


the  system  (36)  takes  the  form: 


Y a = A Ha 


where  the  matrix 


N = K + M 


(38) 


(39) 


has  been  introduced.  The  condition  that  (38)  possesses  other  than  the 
trivial  solution,  = 0 , yields  the  lambda  equation: 


!! 


I!  Y- 


= 0 


(4o) 


for  the  characteristic  frequencies. 

As  the  coordinate  functions  y.  and  OC'  are  not  linearly  dependent, 

0 j 

the  matrix  is  not  singular  and  one  can  presu.lt  iply  both  sides  of  (38) 
by  Y • With: 

\ ~ (41) 

equation  (38)  appears  as: 


where  P = Y 


(42) 


The  associated  characteristic  equation  is» 


- 1 


: s (43) 
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a The  matrix  r is  best  obtained  as  solution 


of  the  matrix  equation: 


where  the  n right  hand  sides  axe  the  n column  vectors  Tlj  • 

weighting  function  method  requires  that  both  sides  of  eq.  (35)  bo  orthogonal 
to  n members  of  a complete  set  of  functions  £ , viz. 

*4  *>,  fi)  = * aJ  f((K,^)  , if.)  + ((/*,«;),  <Pi)}  , m 


With: 


y = Uy)  ; &*(%,<&).  : Kij  •((«,%),  (P{)  , Ks  (KG  < 

V1y»((M,«j),  If,)  ; M = (M.j)  ; N = K*M  , n.j  = (k.j  + M.j) 


this  reduces  to  the  familiar  generalized  eigenvalue  equation: 


y a = ^ Wa 


with  the  associated  frequency  equation  (lambda  equation): 


where 


~ A 


= G . 


Again.,  is  not  sing-alar,  hence  on©  &ay  also  writ©: 

Fa  = A a 
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I with  the  frequency  equations 
8 


- 0 


If  the  first  n members  of  the  complete  set  are  the  functions  1 

x° , x^ , X2,  . . . « , x""1 , 

it  is  said  that  one  requires  equality  between  the  areas  under  the  curves 
representing  the  left  »■  and  right  hand  side  of  (35)  aud  between  their  first 
(n-1)  moments  about  the  y«axis« 

Other  suitable  sets  are  the  coordinate  functions  of  the  Sits  expan- 
sion. If  these  are  orthogonal  a reduction  of  labor  ie  effected.  The  expansion 
functions  which  are  recommended  for  use  are  orthogonal  with  respect  to  the  mass 
distribution  function  J-L  as  a weighting  function.  Their  construction  has 
been  described  in  full  detail  on  pg.  60* 

Belated  to  the  weighting  function  method  is  the  modified  collocation 
method  - outlined  on  pg.  5^»55  - where  it  is  required  that  the  integrals  over 
n adjacent  subintervals  of  both  sides  of  eq.  (35)  be  equal.  This  method  is, 
of  course,  less  cumbersome  than  the  weighting  function  method  itself,  as  the 
integrations  need  not  be  carried  out  over  the  whole  length  of  the  beam,  but 

i 

only  over  n small  portions  of  it.  Numerical  integration  will  have  to  be  used 
in  practically  all  cases,  and  thus  one  integration  process  will  yield  the 
eigenvalue  equations,  if  only  the  cumulative  results  of  this  numerical  inte- 
gration are  registered  at  the  n points.  For,  truly,  the  difference  between 
the  partial  sums  at  and  xt  is  equal  to  the  integral  of  the  function  over 

the  interval  x.  .to  x,  • 

1-1  * 
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f ) On  the  Convergence  of  the  Collocation  Method 

Collocation  methods  have  been  used  in  the  preceding  sections  with- 
out. however,  demonstrating  that  they  lead  to  results  that  approximate  indeed 
the  exact  solution.  Readers  who  are  familiar  with  the  Fredholm  theory  of  in- 
tegral equations  will  know  that  collocation  methods  are  identical  with  the 
classical  treatments  of  linear  integral  equations  "by  Yolterra  and  Fredholm. 

(A  complete  presentation  of  these  methods,  and  in  the  author's  opinion  - the 


clearest  one,  is  given  in  G-.  Kowalewski, 


sn,  Berlin  1930.) 


It  is  instructive  to  outline  the  methods  once  more  end  to  state  the  most  im- 
portant findings  of  the  theory  with  respect  to  application. 

Consider  that  the  following  (non-homogeneoue)  linear  integral  equa- 


tion has  to  he  solved: 


b 

f (>)  ~ J K (*;y)  f (y)  dy  = q(x) 


in  the  interval  as-  x ^ 3.  Let  x,  he  a point  in  the  interval  [a,b  j . 


Then  at  z > the  equation: 


t(*i)  - I K(*i  -,\j)  $(y)dy  = 


is  certainly  satisfied.  Replace  the  integral  by  an  approximating  sum: 


o 

j K(*i  >y)  Uy)Jy  ~ ZL 
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and  1st  the  interval  [a,b]  be  subdivided  into  n equal  intervals,  then  one 
obtains  the  following  n simultaneous  equations; 


f(*0  - ZL  K(Xiiyj)  = g(*;) 

which  are  satisfied  only  approximately  by  the  exact  solutions  of  the  integral 
equation,  x^  and  yj  take  on  all  values  a + k(b-a)/n;  (k  s 1, 2, . . ,n) : 
independently.  With  f;  - ; g;  ~ ; KtJ  = - K (*  i ; )y2"x  , 

one  arrives  at  the  equations:  ^ ■+■  2 K;;  fj  — g«.  or  written  out  in  full: 

(i  4 K"b)  -f,  - K(lfz  v -•  + K,„fn  - g, 


+ Tn 


K^j  f,  + * 11  ’ ^ Kin  ^r» 

Kmf,  •+  Knx-Fx  + •••  + (l  + KM)fn=  S’ 

Cramer's  rule  gives  the  solutions  of  this  system  as: 

■P  - — 

1*  (r  * 1,2,  n) 


where  D is  the  determinant  of  the  coefficients  of  the  unknowns,  and  D is 
obtained  from  D by  replacing  the  r-th  column  in  £ by  { ’ * • * ,Sn  j * 


D = 


K(Z 

0+K,)  1 

iH 

on 

^l.r+i  Km  j 

(i  ♦ • • 

Km 

K*,  - 

• Ks,m 

5* 

Kx}r+ 1 •**  K Hi 

♦ * < . 

• n - 

} r 

Km 

K na 

! 

(l  + 

Km 

K’n.r.. 

3b 

.........  | 

^n.  r*i  K r n | 
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Suppose  tbs  frhavs  been  calculated,  baaed  on  a subdivision  of  the 
interval  [a,b]  Into  n equal  parts,  and,  the  values  of  the  ordinates  are  connect^ 
ed  by  straight  lines,  as  the  following  sketch  illustrates. 

v—- -'T'-'s  The  question  arises  now  as  to  what 

j happens  if  the  subdivision  is  mads 

g ,0 U j£  !■» jn-;  | n 

d i smaller  and  snaller. 


This  necessitates  the  investigation  of  the  so-called  Fredholm  determinant  of 
the  Kernel  K(x,y)  whioh  is  the  limit  of  D as  the  number  n of  subdivisions  of 
[a,b]  tends  to  infinity. 

The  Fredholm  theory  ascertains  the  uniform  convergence  of  the 
approximative  solutions  towards  the  exact  solution  as  the  subdivisions  of  the 
interval  [a*"*5]  become  smaller  and  smaller. 

This  important  result  justifies  the  collocation  method - 
g)  The  Method  of  Station  Functions  (another  method  of  constructing 
coordinate  functions  for  the  RITZ  expansion) 

Qhite  frequently  one  finds  the  suggestion  that  the  so-called  station 
function  method  yields  suitable  coordinate  functions  for  the  RITZ  expansion. 
This  recommendation  will  now  be  examined  critically. 

Recall  that  the  convergence  of  the  RITZ  expansion  is  best,  if  the 
assumed  functions  are  as  cloBe  as  possible  to  the  actual  eigenfunctions  of  the 
problem  under  investigation.  In  section  e)  it  was  suggested  that  the  eigen- 
functions of  a nneighboringB  problem  be  used  to  construct  an  approximation  to 
the  node-loads  of  the  actual  problem,  and  then  to  calculate  the  coordinate 
functions  as  the  static  deflections  under  these  loads.  These  offer  the  great 
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advantage  of  "being  orthogonal  functions,  thereby  reducing  the  numerical  cal- 
culations to  solve  the  matrix  eigenvalue  equation.  But  furthermore  they  are 
already  close  to  the  actual  modes.  The  station  function  method  is  now  dove- 


-LOpea. 


Suppose  the  linear  homogeneous  integral  equation  (Fredholm  equation 


of  the  second  kind)  has  to  be  solved. 


t? 

f(x)  - A J K(x  y)  Uy)  dy  , 


where  the  kernel  K(x,y)  is  given.  As  in  the  Fredholm  theory,  replace  the  in- 
tegral equation  "by  a set  of  n linear  homogeneous  equations.  The  solution  to 
this  set  does  not  satisfy  the  integral  equation  exactly.  Through  the  n points 
an  interpolation  function  can  be  fitted  which  permits  the  determination  of 
functional  values  at  all  points  between  a and  b.  Since  n homogeneous  linear 
equations  admit,  in  general,  n solutions  which  correspond  to  the  n roots  of 
the  characteristic  equation,  one  provides  for  this  fact  by  writing  the  solu- 
tion as  a linear  aggregate  of  n interpolation  functions: 


fc*>  = f(xi)gux)  (52)  I 

Because  f(x)  must  take  on  the  value  f(x^)  for  x = x^,  it  is  apparent  that  the 
functions  g{x)  sure  subject  to  the  requirement: 

= Cji  . [E  = ceu>]  <»> 

The  functions  have  to  satisfy  the  boundary  conditions  and  should*  if  f 

^ 1 

possible,  be  identical  with  the  normal  functions  of  the  problem.  It  is  re-  3 

called  that  the  kernel  function  is  an  "elementary  solution"  of  the  problem  | 

i 

AFRD.n  AQTir  A wrv  ctdt  irn  icrc  Deep  a pru  I 
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from  which,  all  solutions  can  be  built  up  by  means  of  the  superposition  prin- 
ciple. If  h^(z)  is  a load-function,  then 

P 

J K (*;y)  hj (j)  dy 

is  a solution  of  the  problem  which  satisfies  all  boundary  conditions.  Aa  a 
matter  of  fact,  this  integral  represents  the  effect  due  to  the  cause  hj(x). 
Since  every  g^(x)  has  to  satisfy  n conditions,  viz. 

= ~ “ ' ~ (^i* >)  2 0 ) (^C)  w i ~ ~ ^ t 

and  a total  of  n functions  gj_(x)  have  to  b6  constructed,  n^  free  parameters 
must  be  chosen  such  that  all  these  conditions  can  be  fulfilled.  Consequently 
the  g. ’ s must  have  the  form: 


9i(*>  = a<ifK(x;j)hjly)Jy 


I where  the  load  functions  hj(x)  are  given  functions.  They  ought  to  be  approxi- 
mations to  the  mode-loads,  under  which  condition  the  diagonal  elements  a ^ of 

A. 

b (ai 4)  will  be  much  larger  than  the  other  ones. 

J A 

One  is  now  left  with  the  task  of  evaluating  the  matrix  H so  as  to 
I make  g^x^)  s e^.  Let  ^ be  the  vector:  ( 3>  . Si,  , gn  \ > k the 

1 i \ . 

| vector  ...  , with  an  element  £•  - j K (x  / y)  hj  Cy)  dy  ) then  the 


equations  (54)  appear  as: 


(54’) 


If  x takes  on  the  values  x^,  Xg,...,^  this  becomes  a set  of  n simultaneous  I 
equations.  With  the  definitions:  = , they  are  con-  § 

veniently  written  as  ^3  - A K . (55)  I 
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In  view  of  equation  (53)  > 
is  the  reciprocal  of  f'C  ’ 


is  the  unit  matrix  £ , and  thus 


f 

where  d ,j  = (k y)  and  kl}  ~ j K (Xj^y)  (y)  dj  . 

a 

(Note  the  order  of  the  indices!)  Once  the  a^'s  are  known,  one  nay  substi- 
tute eq.  (5L)  into  eq.  (51)  to  obtain: 

fb  - 1 

fCx)—  ^ J K("X ; 'j)  . Z.  . f(*i)  dj  , ^5?) 

a L J 

or  also: 

« b 

fu)  - fooj k(x;>) giCy)dy  . (57. 

J a 

This  equation  has  to  be  satisfied  for  all  x,  especially  for  x s Xj  (j  = 1, 
...,a)  and.  it  possesses  a non-trival  solution,  only  if 

det.  F-~  {■  ■+•  2]f-  feij  j = 0 3 (where  f(xj)  = fj)  (58) 

m 

is  satisfied.  The  roots  of  the  frequency  equation  determine  n sets  of  solu- 
tions £ f ^ | which,  when  substituted  into  (52)  give  a continuous  approxima- 
tion to  the  solutions  of  the  integral  equation. 

Tha  inta,jrat  ions  occuring  in  (58)  will,  in  general,  have  to  be 
carried  out  numerically.  Gaussian  quadrature  yields  the  highest  degree  of 
accuracy  for  a given  number  of  integration  stations.  It  will,  however,  sel- 
dom be  applied,  because  the  values  of  the  function  will  be  given  at  equal 
intervals.  In  any  case,  the  number  of  integration  stations  has  to  be  large 
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enough  bo  render  the  error  of  the  numerical  integration  insignificant  compared 
with  the  approximation  introduced  with  (52). 


It  is  apparent  that  the  station  function  method  is  very  laborious. 

In  the  present  problem  its  application  is  not  Justified,  because  exact  solu- 
tions of  neighboring  problems  can  be  found  readily  and  be  used  in  a RXTZ  expan- 
sion. The  merits  of  the  station  function  method  lie  in  flutter  analysis  of 
aircraft  wings  and  not  in  the  vibration  analysis  of  beams.  There  the  method, 


extended  to  two  independent  variables  (planar  problems)  is  really  useful,  for 
usually  only  the  lowest  frequency  is  of  interest  and  coordinate  functions  to 
be  utilised  in  a R1TZ  expansion  are  not  easy  to  find. 
k.  Limitation  of  the  Analytical  Methods 


In  section  1)  of  this  chapter  five  restrictions  have  been  found  that 
must  be  imposed  upon  the  vibrating  beam  in  order  that  the  present  analysis  be 
applicable  and  it  had  been  stated  that  other  restrictions,  of  dynamical  nature, 
would  be  formulated  later.  To  these  attention  is  now  turned.  Recall  briefly 
the  0 integral  methods'1  that  have  been  outlined,  viz.  the  variational  method 
and  the  integral  equation  method.  The. first  of  these  states  that  the  total 
energy  of  the  vibrating  body  must  be  a minimum  for  every  natural  mode.  The 
second  states  that  the  natural  modes  are  the  characteristic  solutions  of  an 
Integral  equation  which  represents  all  possible  vibrations  of  the  body.  The 
unknown  solutions  are  written  as  linear  aggregates  of  known  vibrations  and 
the  problem  immediately  arises  of  determining  the  degree  of  participation  of 
these  known  vibrations  in  order  to  build  up  the  desired  solutions.  This  leads 
then  tc  an  algebraic  problem  of  solving  a matrix  eigenvalue  equation.  Svery 
solution  vector  corresponds  to  characteristic  solution  of  the  problem.  One 
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I would  be  tamnted  to  believe  that  any  desired  degree  of  accuracy  can  be  obtained! 

j 

if  only  the  Eitz  expansion  (i.s.  the  linear  aggregate)  comprised  a sufficiently^ 

i 

large  number  of  linearly  independent  functions  and  corresponding  part  ideation 
factors*  But  this  is  not  soi  Suppose  that  a Eitz  expansion  contains  twenty 
members  and  that  the  eigenvalue  problem  has  been  solved*  Then  one  may  write 
down  the  expansion  of  the,  say,  tenth  natural  mode  of  the  beam  and  one  knows 
I from  mathematical  reasoning  that  this  ought  to  be  a very  good  approximation 
to  this  mode,  if  the  linear  aggregate  represents  the  solutions  of  a neighbor- 
ing problem.  Let  the  beam  under  investigation  have  a length  of  ten  times  ice 

(average  height.  Under  these  assumptions  a half-wave  cf  the  tenth  mode  has  a 
length  of  about  one  tenth  of  the  span  of  the  beam.  Therefore,  the  nodal 
points  subdivide  the  beam  into  sections  with  a length  to  height  ratio  of 
approximately  one  to  one.  Between  two  adjacent  nodal  points  the  bean  deflects 
as  if  it  were  a portion  of  a continuous  beam  simply  supported  at  the  nodal 
points.  The  simplified  beam  theory  is,  however,  not  applicable  to  such  stubby 
beams,  since  the  stress-strain  distribution  across  the  beam  is  no  longer 
linear.  Shear  deformation  becomes  so  great  that  the  waxping  of  the  cross 
sections  can  no  longer  be  neglected.  Furthermore  the  normal  stresses  perpen- 
dicular to  the  beam  axis  must  be  considered.  This  means  that  a solution  satis- 
fying  the  rigor  of  the  mathematical  theory  of  elasticity  is  called  for.  Un- 

I fortunately  no  such  vibration  solutions  have  been  found  and  it  is  not  possible 
to  compare  the  engineering  solution  with  the  elasticity  solution  directly. 
Bending  problems,  however,  have  b«*er.  solved  in  great  number  and  there  it  was 
found  that  as  long  as  the  beam  has  a length  of  more  than  three  times  its  aver- 

a 

| age  depths,  the  engineering  fosse  las  give  a reasonable  approximation  to  the 
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ij  'behavior  of  the  beam.  The  factor  three  cannot  be  taken  religiously,  for  it 

a 

| depends  on  the  boundary  conditions,  the  geometry  of  the  beam.,  and  the  degree 
| of  accuracy  required.  It  is  reasonable  to  assume  that  the  vibration  analysis 
| obeys  the  same  general  restrictions  as  the  static  deflection  analysis  and, 

I until  more  constraining  requirements  can  be  derived  from  exact  analysis,  it  is 
9 recommended  that  only  those  modes  should  be  calculated  by  the  methods  developed 
I in  this  report  which  yield  half-wave  lengths  that  are  long  in  comparison  to  the 
| average  depth  of  the  beam  ( i. e.  at  least  three  to  five  times  the  depth)* 

The  restriction  just  voiced  should  not  mislead  tha  analyst  to  the 
conclusion  that  is  is  useless  to  carry  many  terms  in  a HIT 2 expansion,  believ- 
ing that  the  higher  terms  represent  the  higher  modes  of  the  beam.  This  is  not 
true,  for,  if  one  had.  an  analytical  expression  in  closed  form,  say  of  the  first 
mode,  and  expanded  this  into  a uniformly  converging  series  and  ueed  the  expan- 
sion instead  of  the  closed  form  in  the  calculations,  then  certainly  all  terms 
of  the  expansion  are  needed  to  represent  the  mode  properly.  The  limitation  of 
the  method  is  of  a different  nature.  It  states  only  that  the  lower  modes  can 
be  calculated  by  the  formulas  given,  because  the  equations  are  valid  only  in 
a domain  of  low  frequencies.  Higher  frequencies  and  modes  do  exist  but  they 
cannot  be  calculated  with  a vibration  theory  that  is  based  upon  a simplified 

I bean  analysis. 

Furthermore,  it  must  be  emphasized  that  the  beams  to  be  Investigated 
have  to  possess  a symmetric  cross  section  and  the  vibrating  has  to  occur  in 

1 

I the  plane  of  symmetry.  Sccentric ( fixed  loads  an  the  be*®  are  thus  to  be 

g s 

S excluded,  as  no  coupling  of  bending  - and  torsional  vibrations  is  permissible.  I 
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Only  recently  B.  Eeilig,  Ingenleur-  Archiv,  vol.  19 * pg*  231-254,  (1951) » has 
developed  a rigorous  solution  for  the  coupled  torsion-tending  vibration  of 
prismatic,  thin-walled,  open-profile  beams,  taking  proper  account  of  the 
warping  stresses.  An  extension  of  this  theory  to  prismatic  bars  of  arbitrary 
cross  section  encounters  great  difficulties  which  arise  from  the  fact  that, 
at  present,  it  is  not  possible  to  calculate  the  stress  distribution  in  & pris- 
matic bar  of  arbitrary  cross  section,  with  arbitrary  boundary  conditions,  due 
to  bending  and  torsion.  Efforts  ought  to  be  made  in  this  direction,  and  the 
above  mentioned  paper  gives  a good  starting  ground. 

The  findings  of  this  report  are  summarized  in  the  following; 

Lemma;  The  methods  presented  in  this  report  are  applicable  to  straight  beams 
with  symmetric  cross  section,  which  vibrate  in  the  plane  of  symmetry.  They 
yield  satisfactory  results  for  modes  whose  half-wave  length  is  long  in  compari- 
son to  the  average  height  of  the  beam. 
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